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ABSTRACT 


Computer  analysis  with  IDL  was  performed  on  polarized  infrared  image  data 
obtained  during  the  MAPTIP  coastal  experiment.  Analysis  of  contrast  enhcuicement 
of  the  oceanographic  ship  Hr.  Ms.  Tydeman  and  a  Lynx  helicopter  using 
polarization  filters  is  presented.  Contrast  of  vertical  and  horizontal  polarized 
images  was  determined  by  two  methods.  Unpolarized  image  contrast  was 
developed  from  polarization  data  and  used  to  calculate  contrast  improvement 
factors  for  vertical  and  horizontal  polarizations.  Background  effects  upon  contrast 
were  explored.  Dependence  of  projected  ship  area  on  contrast  improvement  was 
analyzed.  Vertical  and  horizontal  polarized  images  were  compared  to  verifiy 
preliminary  observations  of  negligible  predominant  polarization  of  man-made 
targets. 

Results  of  this  analysis  confirm  that  vertically  polarized  sea  background 
radiance  can  be  reduced  by  polarization  filtering  to  enhance  ship-to-backround 
contrast.  Man-made  targets  appear  to  present  negligible  polarization.  Sky 
background  immediately  above  the  horizon  also  appears  unpolarized.  Combining 
radiance  from  both  sky  backround  above  the  horizon  and  sea  background  to  form 
an  average  background  reference  radiance  reduces  overall  contrast  improvement. 
Ship-to-background  contrast  should  be  calculated  with  respect  to  the  horizon 
location. 
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I.  INTRODUCTION 


As  the  U.  S.  Navy  moves  to  increase  weapons  systems  capability  in  the  littoral 
environment,  the  ability  to  detect  craft  with  low  or  negligible  radar  cross-section  such 
as  fiberglass  hulled  gunboats  demands  increased  consideration.  Infrared  detection  and 
tracking  systems  enhance  the  combat  capability  of  both  ships  and  aircraft  in  this 
environment.  However,  infrared  signatures  of  interest  often  exist  in  substantially 
cluttered  backgrounds.  Design  of  infrared  systems  capable  of  satisfactory  operation 
despite  clutter  requires  improvement  of  target  to  background  contrast  to  meet  the 
system  detection  threshold  limit.  Successful  clutter  resistant  infrared  systems  could 
fill  the  gap  between  longer  range  radar  detection  and  short  range  daytime  visibility 
detection. 

Current  infrared  detection  and  tracking  systems  use  a  variety  of  filtering 
methods  for  target  detection  and  recognition.  Various  spatial  and  spectral  filtering 
techniques  achieve  some  improvement  in  target  to  background  contrast.  However, 
polarization  filtering  could  provide  significant  improvement  in  target  to  backgroimd 
contrast. 

Polarization  filters  reduce  polarized  apparent  sea  radiance  clutter  in  the  target 
scene  and  improve  target  to  background  contrast.  Specifically,  the  electric  vector  of 
sea  radiance  clutter  appears  perpendicularly  (horizontally)  polarized  within  sim  glint 
regions  and  parallel  (vertically)  polarized  otherwise  (Cooper,  Crittenden,  Milne,  Moss, 
Gregoris,  1992).  Man-made  targets  appear  much  less  polarized  than  sea  backgrounds. 
Polarization  filters  aligned  horizontally  while  viewing  target  scenes  outside  of  sim  glint 
indicate  a  greater  reduction  of  sea  background  radiance  clutter  compared  to  vertically 
aligned  filters  viewing  the  same  target  scene  (Cooper,  Lentz,  Walker,  Chan,  April  94). 
For  infi-ared  system  applications,  the  target  scene  contrast  achieved  with  polarizers 
should  exceed  that  of  an  unpolarized  scene  in  order  to  reduce  the  system  detection 
threshold  value.  This  in  turn  would  increase  minimum  detection  and  recognition 
ranges. 

This  thesis  analyzes  marine  environment  polarization  images  taken  during  the 
MAPTIP  experiment  of  18-31  October  1993  off  the  coast  of  the  Netherlands.  These 
images  provide  the  apparent  radiance  from  target  scenes  involving  ships  and  aircraft 
as  viewed  through  vertical  and  horizontal  filters.  Computer  code  was  generated  to 
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perform  the  analysis,  providing  the  apparent  average  ship  and  background  radiance  in 
vertical  polarization,  horizontal  polarization,  and  impolarized  representation.  Also 
provided  are  vertical,  horizontal,  and  unpolarized  target  to  background  scene  contrasts. 
As  a  result,  contrast  improvement  attained  through  vertical  and  horizontal 
polarizations  compared  with  unpolarized  scenes  is  produced. 

Chapter  II  provides  a  discussion  of  essential  infrared  radiation  fundamental 
concepts  and  equations  necessary  for  more  involved  development  of  infrared  radiation 
theory.  Chapter  III  describes  the  concepts  involved  in  light  polarization  and  the  theory 
of  infrared  polarization  in  the  marine  environment.  Chapter  IV  describes  the  image 
data  utilized  and  also  the  specifics  of  the  analysis  procedure.  Chapter  V  provides  a 
discussion  of  the  results,  conclusions  and  recommendations.  Appendix  A  contains  the 
three  computer  programs  used  and  Appendix  B  provides  the  results  of  the  analysis. 


II.  INFRARED  RADIATION  FUNDAMENTALS 


A.  THEORY  OF  INFRARED  RADIATION 

Infrared  radiation  is  a  part  of  the  total  electro-magnetic  spectrum.  The 
infrared  region  spans  wavelengths  from  0.72  pm  to  1000  pm  as  shown  in  the  diagram 
of  the  electro-magnetic  spectrum  in  Fig.  2.1  (Grum  and  Becherer,  1979).  The  red  end 
of  visible  light  borders  the  infrared  region's  lower  limit  and  the  microwave  end  of  radio 
waves  borders  the  upper  limit.  The  infrared  region  subdivides  further  into  near  (0.72- 
1.5  pm),  middle  (1.5-5. 6  pm),  and  far  (5.6-1000  pm)  infrared  regions  (Hackforth,  1960). 
Any  material  at  a  temperature  above  absolute  zero  emits  radiation.  A  body  above 
absolute  zero  contains  a  finite  value  of  internal  energy  associated  with  the  excitation  of 
atoms  and  molecules.  Excitation  comprises  the  storage  of  energy  in  electronic, 
vibrational,  and/or  rotational  states.  The  levels  of  excitation  in  vibrational  states  occur 
in  multiples  of  some  discrete  value  of  energy.  When  the  atom  or  molecule  returns  to 
the  lower  energy  state,  it  releases  radiation  in  the  form  of  a  photon  at  a  frequency 
proportional  to  the  change  in  energy  between  states, 

A£=Av  2.1 

where,  AE  =  change  in  energy  between  states  (J  or  eV) 
h  =  Plank's  constant  (6.6256  x  10"^  J-s) 
v  =  frequency  of  released  radiation  (Hz) 

The  change  in  energy  associated  with  electronic  states  typically  ranges  from  one  to  ten 
electron  volts  (1.24-0.124  pm),  whereas  vibrational  and  rotational  state  changes  in 
energy  are  t3q)ically  0.1  eV  (12.4  pm)  and  10  ^  eV  (12,424  pm),  respectively.  Therefore, 
molecular  infrared  radiation  primarily  involves  vibrational  changes  of  state  with 
rotational  changes  superimposed.  Fig.  2.2  shows  the  absorption  spectrum  of  diatomic 
molecules  of  HCL  gas.  The  actual  emission  spectrum  contains  rotational  changes  of 
state  superimposed  upon  the  vibrational  changes.  (Cooper  4253  notes,  1994) 
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1994). 


1.  Infrared  Definitions  and  Units 

Several  definitions  and  units  require  introduction  for  an  understanding  of 
infrared  science.  Table  2.1  shows  fundamental  radiometric  quantities  associated  with 
infi'ared  radiation.  Radiant  exitance  (emittance)  or  flux  density,  M,  is  the  radiant  flux 
exiting  an  infinitesimal  area  divided  by  that  infinitesimal  area.  Irradiance  or  radiant 
flux  surface  density,  E,  is  the  radiant  power  incident  upon  a  unit  area  of  surface. 
Radiant  intensity,  I,  is  the  radiant  power  exiting  a  point  source  along  a  given  direction 
within  a  unit  solid  angle.  Radiance  or  radiant  intensity  siirface  density  in  a  given 
direction,  L,  is  the  radiant  power  per  unit  solid  angle  exiting  from  or  incident  on  an 
area  projected  perpendicular  to  the  direction  of  radiant  energy  flow.  (Grum  and 
Becherer,  1979,  Cooper  Class  Notes,  1994;  Wolfe  and  Zissis,  1978,  Bramson,  1968, 
Smith,  1993) 

2.  Blackbody  Radiation 

The  blackbody  is  fundamental  to  the  study  of  infrared  radiation.  Figure  2.3 
shows  t5q)ical  blackbody  radiant  emittance  curves  for  several  temperatures.  The  total 
radiant  emittance  or  area  under  the  curve  and  the  peak  value  of  emittance  depend 
upon  the  blackbody  temperature,  the  peak  value  shifting  to  the  left  at  greater 
temperatures.  The  radiation  emitted  at  a  given  temperature  is  a  maximum  for  a 
blackbody.  As  the  standard  for  radiation  theory,  blackbody  emission  and  absorption 
efficiency,  designated  as  the  emissivity  factor,  equals  one.(Hackforth,  1960) 
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Quantity(S5Tnbol) 

Defining  Equation 

Units 

Radiant  energy  (Q) 

J  (joule) 

Radiant  energy  density  (w) 

cowdQ/dV 

J/m® 

Radiant  power  or  flux  (d>) 

«I>=dQ/dt 

W  (watt) 

Radiant  exitance  (M) 

M=d^>/dA 

W/m^ 

Irradiance  (E) 

E=d«I)/dA 

W/m' 

Radiant  intensity  (I) 

I=d<I>/dco 

W/Sr 

Radiance  (L  or  N) 

L=d^O/dco(dAcos0) 

=dI/(dAcos0) 

W/m=*Sr 

Emissivity  (e) 

e=M/Mbb 

— 

Table  2.1  Radiometric  quantities  (Grum  and  Becherer,  1979). 

CL  Planck's  Law 

During  the  late  19th  century  many  scientists  grappled  with  deriving  the 
equation  fitting  the  experimentally  obtained  radiant  emittance  curves  similar  to  those 
in  Figure  2.3.  Max  Planck  successfully  established  this  equation  and  simultaneously 
ushered  in  the  era  of  quantum  theory.  Planck  theorized  that  oscillating  modes  of 
frequency  co  excite  higher  oscillation  energy  levels  in  discrete  increments  or  quanta. 
This  theory  led  to  Planck's  law  which  follows, 


M,(7)=- 


Inc'^h 

he 

^^-l) 


k^(e 


2.2 


where  Mj^(T)  =  spectral  radiation  emitted  into  a  hemisphere  per  unit  area  per  unit 
wavelength  interval  (Wcm'^cm"^) 

T  =  absolute  temperature  of  the  body  (°K) 

X  =  wavelength  (cm) 
h  =  Plank's  constant  (6.6256  x  10"^  J-s) 

c  =  speed  of  light  in  a  vacuum  (  2.998  x  10^“  cm/s) 

k  =  Boltzmann's  constant  (1.38054  x  10'^^  J-K'^). 
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As  seen  from  this  equation  and  Figure  2.3,  as  the  temperature  increases  the  radiant 
emittance  increases  significantly.  In  fact,  the  total  radiant  exitance  increases  with  the 
fourth  power  of  absolute  temperature  as  shown  by  the  Stefan-Boltzmann  law  described 
below.(Hudson,  1969,  Hackforth,  1960) 


h.  Emiasivity 

No  real  bodies  emit  and  absorb  radiation  ideally  as  the  blackbody. 
Emissivity  depends  on  the  body  constituents,  temperature,  incident  radiation  quantity 
and  wavelength,  and  surface  finish.  The  ratio  of  a  body's  total  radiant  emittance,  M, 
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to  that  of  a  blackbody,  Mbb,  defines  the  body's  emissivity, 


2.3 


Again,  blackbody  emissivity  equals  one.  Some  materials  approximate  blackbody 
behavior.  Many  others,  called  gray  bodies,  display  emissivity  factors  less  than  one,  but 
are  independent  of  wavelength  as  is  a  blackbody.  Still  other  materials,  called  selective 
emitters,  emit  and  absorb  radiation  variably  dependent  upon  wavelength.  The 
following  equation  accounts  for  the  possible  spectral  dependence  of  emissivity,  e(A,), 


€  = 


jeik)M^dk 

0 


2.4 


Figure  2.4  displays  radiant  emittance  curves  for  a  gray  body  and  selective  emitter 
compared  with  a  blackbody  curve.(Hudson,  1969) 
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Figure  2.4  Radiant  emittance  for  selective,  graybody,  and 
blackbody  emission  spectra  (Hudson,  1969). 
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c.  Stefun-BoUzmann  Lmw 

The  integration  of  Planck's  Law  over  all  wavelengths  produces  the  total 
radiant  exitance  emitted  from  a  blackbody  as  a  function  of  temperature  as  shown 
below, 


M(J)=jMf(iX)dX=aJ*[Wlcm^]  2.5 

0 

where,  M(T)  =  total  radiant  emittance  at  temperatiu-e  T  (Wcm'^) 
a  =  Stefan-Boltzmann  constant  (5.670  x  10'^^  Wcm'^K"^). 

For  a  graybody,  this  function  becomes, 

MiT)=eaT*.  2.6 

Quite  often,  spectral  filters  used  with  infrared  systems  permit  only  a  specified 
wavelength  interval  of  irradiance  to  reach  the  detector.  Therefore,  effective  radiant 
emittance  then  results  from  integrating  M(T)  over  only  the  desired  wavelength 
interval.(Hudson,  1969,  Hackforth,  1960) 

d.  Wien's  Displacement  Law 

Differentiating  the  radiant  emittance  of  Equation  2.2  and  setting 
the  result  equal  to  zero  yields  the  wavelength  of  maximum  emittance.  This  is  known 
as  Wien's  displacement  law, 

A.„r=K  2.7 

where,  =  wavelength  of  maximum  emittance  (cm  or  pm) 

T  =  absolute  temperature  (°K) 

K  =  constant  =  0.2897  cmK  for  in  cm  =  2897  pmK  for  in  pm. 

Since  the  wavelength  of  maximum  emittance  times  the  absolute  temperature  equals  a 
constant,  then  A„,  displaces  to  shorter  wavelengths  as  temperature  increases  as  shown 
in  Figure  2.3  by  the  dashed  line.(Hudson,  1969,  Hackforth,  1960) 

e.  Kirchoffs  Law 

A  blackbody  absorbs  all  incident  radiation  and  emits  an  equal  amovmt  of 
radiation.  Applying  conservation  of  energy  to  a  non-blackbody  yields, 
where  Uj  =  incident  radiant  energy 
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u,=u  +u  +u 

i  P  O  T 


2.8 


Up  =  reflected  radiant  energy 
U„  =  absorbed  radiant  energy 
=  transmitted  radiant  energy. 

Normalizing  with  respect  to  Uj  produces 

l=p+a+T  2.9 


where,  p  =  Up/U;  =  reflectivity 

a  =  UyUj  =  absorptivity 
X  =  U/Uj  =  transmissivity. 

Kirchoff  showed  that  the  ratio  of  radiant  emittance,  M,  to  absorptivity  for  a  body  at  a 
specified  temperature  equals  a  constant  value  of  the  radiant  emittance  of  a  blackbody, 
MgB,  at  the  same  temperature  or 

2.10 

a. 

Substituting  in  the  Stefan-Boltzmann  equation  produces 

2.11 

a 

a=e.  2.12 

So,  the  emissivity  of  a  body  equals  its  absorptivity.  For  opaque  bodies,  x  =  0,  thus 

e=l-p  2.13 

which  means  good  emitters  reflect  poorly.(Hudson,  1969) 
f.  Lombert's  Lcav 

Since  a  blackbody  is  a  diffuse  radiator,  it  has  the  same  radiance  in  all 
directions.  However,  the  radiant  intensity  of  an  element  of  the  blackbody  siuface 
varies  with  the  angle  from  the  normal  to  the  surface.  Lambert's  Law  states  this  effect 
as. 
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dl^LdAcasB 


2.14 


where,  dl  =  element  of  radiant  intensity 
L  =  radiance 

dA  =  siu-face  area  element  of  body 
0  =  angle  from  normal  to  observation  direction. 

Figure  2.5  shows  the  geometry  involved  in  Lambert's  Law.(Jamieson,  Plass,  McFee, 
Grube,  and  Richards,  1963,  Hackforth,  1960).  The  n  steradians  in  Figure  2.5  denotes 
the  quadrant  of  space  to  the  left  of  the  source  area  and  above  the  horizontal  line  that 
passes  through  the  center  of  the  source. 


B.  INFRARED  RADIATION  SOURCES 

As  stated  before,  any  body  above  absolute  zero  emits  radiation.  Types  of 
infrared  radiation  sources  include  artificial,  natural,  and  background.  Figure  2.6 
depicts  the  many  t3q)es  of  infrared  radiation  sources.  Artificial  somces  consist  of 
primarily  man-made  materials  such  as  lamps,  filaments,  and  laboratory  blackbody 
sources.  Trees,  water,  mud,  vehicles,  and  people  makeup  terrestrial  natural  sources. 
Atmospheric  gas  molecules,  dust,  and  water  droplets  are  also  terrestrial  sources. 
Terrestrial  sources  emit,  reflect,  absorb  and  reradiate  radiation  characteristic  to  each 
particular  material.  Celestial  natural  sources  include  the  sim,  moon,  planets,  and 
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stars.  Terrestrial  and  celestial  sources  serve  as  background  source  radiation  if  super¬ 
imposed  on  the  primary  source  radiation.(Hackforth,  1960) 


Background  radiation  also  differs  with  respect  to  day  or  night.  Figure  2.7 
depicts  the  relative  radiant  energy  of  clear  day  and  night  skies.  During  the  day, 
reflection  and  scattering  of  sunlight  by  atmospheric  particles,  water  droplets,  dust 
particles,  and  clouds  as  well  as  emission  from  atmospheric  constituents  at  ambient 
temperatures  contribute  to  backgroimd  radiation.  Da3d;ime  backgroimd  radiation 
increases  in  the  direction  of  the  sun  and  horizon  and  remains  higher  than  nighttime 
levels  in  all  directions.  The  angle  of  elevation  and  atmospheric  temperature  affect  both 
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Figure  2.7  Spectral  energy  distribution  of  background  radiation  from  clear  day  and  night 
skies  (Hackforth,  1960). 

daytime  and  nighttime  radiation  intensity.  At  night,  the  high  frequency  radiation 
scattered  by  small  particles  and  molecules  vanishes.  Background  radiation  at  night  is 
emitted  from  atmospheric  molecules  and  particles  with  the  moon  and  stars  providing 
additional  low-level  radiation.  (Hackforth,  1960) 

Discrimination  of  a  primary  source  against  its  background  requires  a  finite 
primary  source  temperature  different  from  its  background,  or  a  difference  in  source 
and  background  emissivities.  The  magnitude  of  the  required  temperature  difference 
depends  upon  both  the  detection  system  capabilities  and  the  ambient  photon  noise 
environment.  The  background  signal  generates  noise  for  the  detector.  Systems  use 
various  methods  such  as  spatial,  spectral,  or  electronic  filtering  to  reduce  the 
background  signal  to  enhance  primary  source  or  target  detection.  (Hackforth,  1960) 

C.  ATMOSPHERIC  PROPAGATION  OF  INFRARED  RADIATION 

Whatever  the  source,  the  detectable  radiation  received  at  a  detector  differs  from 
that  radiation  emitted  at  the  source.  The  atmosphere  or  some  other  intermediate 
medium  through  which  the  radiation  passes  attenuates  the  radiation  through 
scattering  and  absorption.  These  effects  occvu-  independently  from  the  source,  but  vary 
according  to  the  wavelength  of  radiation.  Current  meteorological  conditions  play  an 
important  role  in  radiation  propagation.  As  the  atmosphere  changes  so  does  its 
propagation  characteristics.  Cloud  type,  base,  and  ceiling,  fog,  suspended  particles, 
and  molecular  constituents  of  the  atmosphere  all  impact  absorption  and  scattering 
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effects  of  the  atmosphere  on  infrared  radiation.  Beer's  Law  quantifies  the  total 
attenuation  by  the  atmosphere  as 


I{z)=me-'^  2 

where,  1(0)  =  initial  intensity  of  radiation 

I(z)  =  intensity  of  radiation  after  propagation  distance  z  in  the 

atmosphere 

p  =  extinction  coefficient  of  atmosphere 
z  =  propagation  distance. 

Therefore,  the  transmittance,  t  results  from, 

x=M.=e-^.  2. 

m 

The  atmospheric  absorption  (p^)  and  scattering  (p^)  coefficients  sum  to  produce  the 
total  atmospheric  extinction  coefficient  (p) 


where. 


m  a 


u,=o_+o„  2.19 

and  K„,  =  molecular  absorption  coefficient 
Ka  =  aerosol  absorption  coefficient 
=  molecular  scattering  coefficient 
Og  =  aerosol  scattering  coefficient. 

The  magnitudes  of  each  coefficient  depend  upon  the  atmospheric  composition  and 
constituent  density  and  distribution.  Aerosol  scattering  dominates  infrared  radiation 
extinction  below  1.2  pm.  Above  1.2  pm,  molecular  absorption  dominates.  In  addition, 
atmospheric  transmittance  depends  significantly  on  wavelength.(Cooper  class  notes, 
1994)  Figure  2.8  shows  transmittance  as  affected  by  infrared  absorption  by  the 
predominant  atmospheric  molecules  over  a  horizontal  path  of  6000  meters  at  sea  level 
with  17  mm  of  precipitable  water  (Hudson,  1960). 
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1.  Absorption  Effects 

Absorption  occurs  when  infrared  radiation  interferes  with  a  free  molecule,  atom 
or  a  liquid  or  solid  at  about  the  frequency  corresponding  to  an  energy  level  separation, 
and  the  electric  field  of  infrared  radiation  excites  that  material  to  the  higher  energy 
level.  The  radiation  disapears  as  a  result  (Hackforth,  1960).  Lorentzian  curves  depict 
the  frequency  lineshape  of  absorption  which  follows  the  classical  resonance 
phenomenon.  Inhomogeneous  broadening  results  in  frequency  lineshape  that  are 
Gaussian.  Also,  collisional,  Doppler,  or  pressiu-e  broadening  increases  the  possible 
range  of  frequencies  for  absorption  (Cooper  class  notes,  1994).  Carbon  dioxide  and 
water  dominate  infrzured  molecular  absorption  in  the  atmosphere  (Hudson,  1969). 

To  a  much  lesser  extent,  absorption  and  scattering  by  aerosols  affects  infrared 
propagation.  Particles  suspended  in  the  atmosphere  such  as  dust  and  volcanic  ash 
generate  aerosols.  Their  distributions  and  number  densities  vary  considerably  from 
one  environment  to  another.  Several  aerosol  models  available  predict  aerosol 
composition  of  the  atmosphere  and  effects  on  propagation.(Cooper  class  notes,  1994) 

2.  Scattering  Effects 

Scattering  results  from  the  absorbing  molecule,  atoms,  or  aerosol  subsequently 
reradiating  but  over  a  4n  solid  angle,  thus  reducing  the  directivity  of  the  original 
radiation.  Where  the  incident  radiation  wavelength  greatly  exceeds  the  size  of  the 
molecule  or  aerosol  particle,  the  following  reradiation  approximates  isotopic  radiation. 
As  particle  size  increases  compared  to  the  incident  radiation  wavelength,  the  radiation 
increasingly  tends  to  be  scattered  into  the  forward  hemishpere.  Rayleigh  scattering 
theory  effectively  describes  small  particle  scattering,  while  Mie  theory  covers  the  entire 
range  of  particle  sizes.(Cooper  class  notes,  1994) 
a.  Rayleigh  Scattering 

Rayleigh  scattering  theory  describes  scattering  for  small  particles  of 
radius  r,  typically  molecules,  where  2kt/'K«1.  Other  assumptions  for  Rayleigh 
scattering  include  nonionized  isotropic  particles,  index  of  refraction  approximately  one, 
and  incident  radiation  whose  frequency  is  different  from  particle  resonant  frequencies. 
Since  diatomic  nitrogen  and  oxygen  constitute  most  of  the  atmosphere,  the  isotropic 
assiunption  requires  a  molecular  anisotropy  correction  factor.  The  equation  governing 
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Rayleigh  scattering  follows, 


-  =K'*— (1  +cos^Q)  2.20 

I' 

where,  I  =  intensity  of  scattered  radiation 
r  =  intensity  of  incident  radiation 
K  =  characteristic  of  medium 
V  =  particle  volume 

r  =  distance  to  particle  from  observation  point 
A,  =  wavelength  of  radiation 
6  =  scattering  angle. 

S5mchronously  induced  dipole  moments  from  the  incident  radiation  produce  Rayleigh 
scattering.  (Smith,  1993)  The  scattered  radiation  frequency  equals  the  incident 
radiation  frequency.  Thus  Rayleigh  scattering  is  elastic  (Coulson,  1988). 
b.  Mie  Scattering 

Mie  scattering  applies  to  both  large  and  small  particles  compared  to 
incident  wavelength.  This  theory  assumes  isolated  spherical  particles  as  the  scattering 
medium.  Although  most  particles  are  non-spherical,  large  numbers  of  randomly 
orientated  non-spherical  particles  approximate  spherical  particle  distributions  except 
for  determining  properties  of  backscatter  and  polarization.  (Smith,  1993)  As  with 
Rayleigh  scattering,  Mie  scattering  describes  the  excitation  of  particles  by  incident 
radiation  and  the  reradiation  of  that  energy  into  a  An  steradisui  solid  angle  (Coulson, 
1988).  The  Mie  scattering  cross  section  per  unit  volume  follows, 

so 

0 

where,  (3^^  =  scattering  cross  section  per  unit  volume 

Csca(r)  =  scattering  cross  section 

n(r)  =  particle  size  distribution  function  (#  particles/unit  volume 
between  r  and  r+dr  (Smith,  1993) 
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III.  INFRARED  POLARIZATION  IN  THE  MARINE  ENVIRONMENT 


A.  POLARIZATION  OF  LIGHT 

Much  experimentation  by  Fresnel  and  Arago  led  to  the  conclusion  that  a  light 
wave  consists  of  two  components  perpendicular  to  the  direction  of  propagation.  The 
components  for  a  wave  propagating  in  the  z-direction  may  be  represented  by 

£/z,0=Fa,cos(or+6,)  3.1 


Ey(z,t)=E^cos(<j»t+&^  3.2 

where  and  E(,y  denote  the  component  electric  field  amplitudes,  5,  and  5y  are  the 
respective  phases,  (o  is  frequency,  and  t  is  time.  E„(z,t)  and  Ey(z,t)  produce  a  vector  in 
the  plane  perpendicular  to  propagation  that  generates  a  curve  in  that  plane  known  as 
the  polarization  ellipse  described  by 

p  2  p  2  EE 

— 3L-2— — ^cos6=sin^6  3.3 

V  V 

where  8  =  8y  +  8*.  This  ellipse  is  shown  in  Fig.  3.1.  Degenerate  forms  of  the 
polarization  ellipse  include  circular  and  linear  polarizations.(Collett,  1993) 


Figure  3.1  Polarization  Ellipse  (Collett,  1993). 
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The  polarization  ellipse  is  applicable  to  completely  polarized  light  only,  and  not 
useful  for  describing  partially  or  unpolarized  light.  Sir  Gleorge  Gabriel  Stokes  resolved 
this  problem  by  describing  light  with  observable  intensities.  He  determined  that  light 
of  any  polarization  can  by  described  by  the  four  Stokes  parameters  given  by 


3.4 

3.5 

Sj  “2£’q^qj,cos6 

3.6 

3.7 

where 

II 

3.8 

(Collett,  1993) 

Then  the  beam  intensity  is  given  in  terms  of  Stokes  parameters 

as. 

7(0, (|)) =— [Sq +5jCOs(20)  +52cos((|>)sin(20)  +53sin(<|>)sin(20)] 

3.9 

where  <|)  is  an  induced  phase  angle  (Collett,  1993). 

Stokes  parameters  of  two  independent  light  beams  may  be  added  to  obtain 
Stokes  parameters  for  the  resultant  beam.  Thus  partially  polarized  light  may  be 
represented  by  the  summation  of  unpolarized  light  with  completely  polarized  light. 
Both  xmpolarized  light  and  partially  p  iarized  light  may  be  represented  by  summing 
two  oppositely  polarized  beams.  For  impolarized  light,  the  intensities  of  the  two  beams 
must  be  equal.  Thus  a  partially  polarized  target  scene  in  a  marine  environment  may 
be  represented  by  two  beams  of  vertical  and  horizontal  polarization.(Collett,  1993) 

B.  REFLECTION  AND  REFRACTION 

Reflection  and  refraction  of  light  is  described  by  the  Fresnel  equations  and 
Snell's  law.  The  incident,  reflected,  and  transmitted  waves  are  shown  in  Figure  3.2. 
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Serving  as  a  starting  point  are  Maxwell's  equations  of  electro-magnetism  in  a  source- 
free  medium, 

3.10 

dt 


vxH=e'—  3.11 

dt 

v»£=0  3.12 

v^=0  3.13 

where  E  and  H  are  the  electric  and  magnetic  field  intensities,  respectively.  The 
complex  permitivity  and  complex  permeability  are  denoted  by  e*  and  p*,  respectively. 
For  plane  waves  of  the  form. 


3.15 

Y=ya)(|i*e*)''*=o+7P  3.16 

incident  at  a  bovmdary  between  two  media  of  different  complex  refractive  indices  nj* 
and  n2*  the  angle  of  reflection  equals  the  angle  of  incidence.  The  angle  of  transmission 
0t  as  a  fimction  of  incident  angle  0;  is  given  by  Snell's  law 

ni*sin(0.)=/i2*sm(0,).  3.17 
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For  nonmagnetic  media  the  Fresnel  equations  can  be  written 


sin(9^-ep 
^  sm(0,+0,) 


tan(0,-0,) 
'  tan(0,+0,) 


2sin(0,)cos(0i) 

t= - - - - 

sm(0j+02) 


2sin(0,)cos(0,) 
sm(0j+0,)cos(0, -0,). 


3.18 


3.19 


3.20 


3.21 


The  Brewster  or  polarizing  angle  occirrs  where 


3.22 


and 


0,+0,=9O. 


3.23 


The  reflection  coefficients  then  become 


r,=0 

3.24 

r,^0. 

3.25 

Thus  the  reflected  light  is  completely  polarized  normal  to  the  plane  of 
incidenceXSandus,  1965) 
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C.  SEA  SURFACE  INFRARED  POLARIZATION 


Background  infrared  radiation  from  the  sea  surface  is  composed  of  reflected 
radiation  from  the  sky,  sim,  moon,  stars,  etc.  Another  contribution  to  background 
radiation  comes  from  thermally  emitted  radiation  from  the  water,  and  molecular  and 
particle  volume  scattering  (Sidran,  1981).  Sea  siuface  reflection  is  governed  by  the 
Fresnel  equations  with  polarization  of  reflected  radiation  greater  near  the  Brewster 
angle.  Polarization  of  emitted  radiation  is  somewhat  more  involved. 

1.  Sea  Surface  Reflectance  and  Emissivity 

Polarized  light  emanating  from  a  siuface  may  be  due  to  emitted  radiation  as 
well  as  reflected  radiation.  Arago  showed  that  light  from  an  incandescent  solid  or 
liquid  is  partially  polarized  even  for  non-polished  surfaces.  Arago's  work  was 
qualitative.  Millikan  later  described  Arago's  work  in  quantitative  terms.(Sandus, 
1965) 


23 


Millikan  studied  the  polarization  characteristics  of  uranium  glass  fluorescence. 
He  observed  that  uranium  glass  fluorescence  emission  reached  a  maximum 
polarization  at  near  grazing  angles.  Furthermore,  the  emitted  light  was  polarized 
parallel  to  the  plane  of  emission.  As  described  above,  at  least  partial  polarization 
should  be  observed  for  reflected  light.  However,  this  light  would  be  polarized  normal 
to  the  plane  of  incidence.  Therefore,  the  polarized  light  Millikan  observed  from  the 
fluorescence  of  uranium  glass  could  only  have  come  from  emitted  light  imdergoing 
refraction.(Sandus  1965) 

Extending  the  understanding  of  emission  polarization,  Basener  and  McCoyd 
developed  a  model  to  predict  the  percent  polarization  of  emitted  infrared  radiation 
from  the  sea  surface  as  a  function  of  wind  speed  and  viewing  angle.  This  model  is 
based  upon  a  sea  surface  of  tiny  plane  facets  of  varying  orientations.  The  facet 
orientation  is  statistically  distributed  depending  on  wind  speed.(Basener  and  McCoyd, 
1967) 


Figure  3.3  Reflection  geometry  for  a  facet  of  a  wind  roughened  sea  surface  (Gregoris,  Yu, 
Cooper,  Milne,  1992). 


24 


Given  that  the  sea  surface  emits  polarized  radiation,  Basener  and  McCoyd 
investigated  the  affects  on  emissivity  of  the  variable  nature  of  the  sea  surface.  Figure 
3.3  shows  the  reflection  geometry  for  a  facet  of  a  wind  roughened  sea  surface.  Each 
facet  is  assumed  flat.  For  a  flat  blackbody  surface  z=f(x,y)  the  directional  emissivity 
e,(9,(l))  is  related  to  the  directional  reflectance  R(0,(j))  through 

e,(0,<|))=l-/?/e,<|))  3.26 

where  r  is  polarization  mode.  Thus,  the  radiant  flux  emitted  per  unit  solid  angle  per 
unit  area  per  unit  wavelength  interval  is 

F,(0,<l.)=|€/e,<l))cos(e)B(X.7)  3.27 

where  B(A,,T)  is  the  Planck  blackbody  fimction.  The  reflectance  is  found  from  the 
Fresnel  equations  above.  For  complex  indices  of  refraction,  the  emissivity  is  given  by 


where  z  ^  is  a  complex  index  of  refraction.(Basener  and  McCoyd,  1967) 

The  above  equation  is  for  a  flat  surface  and  the  sea  surface  is  not  flat. 
However,  the  sea  surface  can  be  approximated  by  many  facets  with  distributed 
orientations.  For  the  facets  to  be  considered  flat  then  they  must  be  large  compared  to 
the  wavelength  of  light.  Sea  surface  waves  easily  meet  this  constraint.  Referring 
again  to  Fig  3.3,  the  vector  O  points  in  the  direction  of  the  observer,  S  is  the  sun 
direction,  and  the  vector  N  is  the  facet  normal.  Therefore,  the  emissivity  depends  on 
the  angle  between  S  and  N.  The  radiant  flux  emitted  per  unit  solid  angle  per  unit 
wavelength  interval  from  a  unit  horizontal  area  is  given  by 

F,(0.4>)=|b(A,7)  f e;(P.Y)sec(Pl/i:p,Y)cos(i|r>ia)  3.29 

where  e'(P,Y)  is  the  emissivity  referenced  to  the  main  plane  of  vision  defined  by  S  and 
k,  where  k  is  the  unit  vector  in  the  z-direction,  P  is  the  zenith  angle  of  the  facet 
normal,  and  y  is  the  azimuth  of  the  facet  normal.  The  facet  distribution  function  is 
denoted  by  f(P,y)  and  i|/  is  the  angle  between  S  and  N.  Cox  and  Munk  (Basener  and 
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McCoyd,  1967)  have  determined  experimentally  that  the  distribution  of  slopes  given  by 
f(P,y)  as  a  fimction  of  wind  speed  is  given  by 


/p,Y)do)=p(z,(P,Y)^(P,Y)) — \ — du>  3.30 

^  cos^CP) 

where 

z^=-tan(P)sin(Y)  3.31 

Zy=-tan(P)cos{Y)  3.32 

and  denote  the  surface  slope  where  y  points  upwind  and  p(z,y)  gives  the  fraction  of  a 
horizontal  unit  area  of  siu-face  whose  slopes  are  within  z^  ±  1/2  5z,  and  Zy±l/25Zy 
(Basener  and  McCoyd,  1967).  The  facet  slope  distribution  for  crosswind  and  upwind- 
downwind  situations  are  shown  in  Figures  3.4  and  3.5,  respectively.  The  horizontal 
gyi-q  shows  the  value  of  beta,  the  surface  normal  angle  from  the  vertical  direction 
(Sidran,  1981). 

The  above  formulation  ignores  affects  due  to  shadowing.  Correction  factor 
equations  were  developed  by  Wagner  and  also  by  Saunders  (Sidran,  1981  and  Gregoris, 
Yu,  Cooper,  and  Milne,  1992).  Omitting  spectral  dependencies  and  including  Savmders' 
correction  factor  S(9),  the  average  apparent  radiance  from  the  sea  is 


2cos(0^) 


3.33 


tan(e„)cos(e,,Y„.9j,(|>^9„.Y,^ 

cos*(9„) 


.(l-e^^9>^9p)iV„^;(9^<|)^9,)l. 
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SLOPE  PROBABILITY  (ARBITRARY  UMTS) 


Figure  3.4  Facet  slope  distribution  -  crosswind  case  (Sidran,  1981). 
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Figure  3.5  Facet  slope  distribution  -  upwind/downwind  case  (Sidran 


Then  the  effective  polarization  is  represented  by 


3.34 


where  Np^y,  is  the  unpolarized  path  radiance  which  can  be  obtained  with  LOWTRAN 
atmospheric  modeling.  Fig  3.6  shows  the  theoretical  values  of  polarization  for 
MWIR(4-5.5vun)  and  LWIR(8-1 1.5pm)  in  and  near  sim  glint  regions.  Note  the  negative 
values  of  polarization  outside  of  the  sun  glint  regions.  Here  the  sea  emission  is  greater 
than  the  reflected  sky  contribution.  The  negative  values  indicate  parallel  or  vertical 
polarization.  The  effect  is  greater  in  the  LWIR  since  the  solar  radiance  is 
predominantly  in  the  MWIR.  (Gregoris,  Yu,  Cooper,  Milne,  1992) 

Also,  if  Q  represents  the  percent  polarization  of  sea  surface  emitted  radiation, 

then 

(j= 100^1^  3.35 

e.+e,. 


Figure  3.7  shows  emissivities  and  percent  polarization  as  a  function  of  viewing  angle  0. 
Polarization  is  independent  of  surface  roughness  for  theta  less  than  approximately  25 
degrees.  Beyond  this  value  of  theta,  the  polarization  increases  with  viewing  angle,  but 
is  partially  reduced  with  increasing  roughness.  In  this  case  an  increase  in  polarization 
again  indicates  parallel  polarization.  Figure  3.8  shows  observed  data  by  Hale  and 
Querry  for  emissivities  and  percent  polarization  for  both  reflected  and  emitted 
radiation  for  all  wavelengths  below  10  microns.  The  observed  and  theoretical  results 
are  comparable.(Sidran,  1981) 
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Figure  3.7  Emissivities  and  percent  polarization  as  a  function  of  viewing  angle 
(Sidran,  1991). 
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Figure  3.8  Emissivities  and  percent  polarization  for  reflected  and  emitted 
radiation  for  all  wavelengths  below  10  microns  (Sidran,  1991). 


PERCENT  POLARIZATION 


2.  Infrared  Polarization  of  Ship  Signatures  and  Background  Contrast 

As  stated  before,  the  infrared  radiance  from  the  sea  surface  is  composed  of  two 
components.  These  are  the  thermally  emitted  and  reflected  radiation  given  by 


3.37 

where  each  component  depends  on  the  reflectivity  and  emissivity  of  the  sea  surface 
averaged  over  the  distribution  of  wave  facet  orientations.  The  emissivities  for  parallel 
and  perpendicular  polarizations  are 

e^6^cos^(a)+ej;sin^(a) 


e^€*sin2(a)+e^cos2(a). 

In  terms  of  average  apparent  radiance  the  degree  of  polarization  is  given  by 

<N>-<Nt.> 

POL=- 


3.39 


<N>+<N^> 


3.40 


and  the  radiance  difference  (absolute  contrast)  for  any  polarization  can  be  written  as 

Thus,  the  polarization  improvement  factor  for  contrast  for  horizontal  filters  over 
vertical  filters  is 


AAT-AiV. 

V _ « 

AN^+ANj 


3.42 


(Cooper,  Lentz,  Walker,  Chan,  1994). 


Contrast  relative  to  the  background  is  most  commonly  defined  as  the  difference 
between  target  radiance  and  background  radiance  divided  by  the  backgrotmd  radiance, 

343 

^btgd- 

Contrast  improvement  then  is  defined  as  the  ratio  between  polarized  contrast  and 
unpolarized  contrast, 

Contrast  between  an  impolarized  target  against  a  sea  background  outside  of 
sun  glint  regions  should  be  improved  with  the  use  of  polarization  filters.  Man-made 
targets  do  not  indicate  any  predominant  polarization  (Cooper,  Lentz,  Walker,  Chan, 
1994).  Therefore,  horizontal  polarization  filters  would  reduce  the  predominantly 
vertically  polarized  sea  background  radiance.  Although  each  filter  will  reduce  the 
unpolarized  target  radiance  by  a  factor  of  approximately  two,  the  reduction  in  sea 
background  radiance  should  be  greater.  This  will  increase  the  ntunerator  of  Equation 
3.43  for  horizontal  polarization  while  reducing  the  denominator.  Thus,  an  overall 
improvement  in  contrast  would  be  realized. 
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IV.  DATA  ANALYSIS 


A.  MAPTIP  EXPERIMENT 

The  polarization  image  data  utilized  in  this  thesis  was  obtained  during  the 
Marine  Aerosol  Properties  and  Thermal  Imager  Performance  Trial  (MAPTIP)  held 
October  11  through  November  5,  1993.  The  purpose  of  the  experiment  was  to  obtain 
atmospheric  propagation  data  below  10  meters  elevation  in  the  3-5pm  and  8- 12pm 
bands.  Acciu-ate  atmospheric  prediction  codes  are  necessary  for  properly  analyzing  the 
performance  of  electro-optical  detection  and  tracking  systems.  Ciurently  the  Navy 
Atmospheric  Model  (NAM)  used  in  LOWTRAN7  effectively  predicts  infrared 
propagation  in  the  marine  boimdary  layer  above  10  meters.  The  Navy  Oceanic 
Vertical  Atmospheric  Model  (NOVAM)  predicts  d5mamic  marine  aerosol  composition  for 
specific  marine  environments  within  the  marine  boimdary  layer.  MAPTIP  is  expected 
to  provide  coastal  atmospheric  propagation  information  to  extend  NOVAM  to  the 
littoral  environment  (Jensen  and  de  Leeuw,  1993).  Image  data  used  in  this  thesis  was 
obtained  during  the  experiment  and  is  of  the  Hr.  Ms.  Tydemann  (TDM),  an 
oceanographic  ship,  and  of  a  L5tix  helicopter. 

B.  DATA  COLLECTION  EQUIPMENT 

Image  data  was  acquired  during  MAPTIP  using  the  NPS  AGA  780 
Thermovision  Thermal  Imaging  system.  Images  were  stored  using  AGEMA  CATS  data 
acquisition  program. 

1.  AGA  780  Thermovision  Thermal  Imaging  System 

The  version  of  the  AGA  780  Thermovision  Thermal  Imaging  System  used 
during  MAPTIP  is  a  dual  band  serial  scanning  infrared  imager.  The  AGA  780 
operates  two  scanning  systems  side-by-side.  One  system  covers  the  3-5.6  pm  and  the 
other  covers  8-14  pm  wavelength  bands.  Each  utilizes  a  7°x7°  field  of  view  lens  that 
focuses  incoming  infrared  radiation  through  a  dual  rotating  prism  system  onto  their 
respective  detectors.  The  short  wave  lens  is  silicon,  and  its  detector  is  an  indium 
antimonide  (InSb)  photovoltaic  detector.  The  long  wave  lens  is  germanium  and  its 
detector  is  a  mercury-cadmium-telluride  (HgCdTe)  photovoltaic  detector.  The  detectors 
are  cooled  to  77K  by  dewars  filled  with  liquid  nitrogen.  The  radiation  passes  through 
a  vertically  mounted  prism  rotating  at  180  rpm  followed  by  a  horizontally  mounted 
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prism  rotating  at  18000  rpm.  The  rotating  prisms  are  synchronized  electronically  to 
the  monitor's  vertical  and  horizontal  scans.  This  S5Tichronization  of  the  prisms  serves 
to  generate  four  fields  of  100  horizontal  scan  lines  each.  Only  70  lines  in  each  field  are 
used  for  producing  an  image;  the  four  fields  are  interlaced  to  produce  an  image  of  280 
lines.  The  scanner  operates  at  25  fields  per  second.  Thus,  6.25  images  per  second  are 
generated.  The  detectors  produce  a  current  signal  called  "thermal  value"  proportional 
to  the  intensity  of  radiation  received.  This  signal  is  amplified  and  converted  into  a 
video  signal  for  monitor  display.(AGA  Operating  Manual,  1980) 

2.  Thermal  Imager  Calibration  Technique 

The  Stefan-Boltzmann  law  predicts  a  bolometric  intensity  proportional  to  the 
fourth  power  of  temperature.  Although  the  intensity  of  radiation  received  at  the 
detector  is  presumed  to  be  proportional  to  the  system's  "thermal  values,"  this  signal  is 
a  non-linear  function  of  temperature.  The  AGA  780  requires  calibration  with  a 
laboratory  blackbody  source.  The  system  output  of  "thermal  value"  is  recorded  against 
the  blackbody  temperature  through  a  range  of  temperatures.  The  calibration  curve  is 
generated  through  a  least  squares  fit  method  and  is  of  the  form 


4.1 


where  A,  B,  C  are  constants  specific  to  the  combination  of  scanner,  lens,  and  filters 
used  during  the  calibration.(AGA  Operating  Manual,  1980)  The  OC  is  for  an  offset 
correction.  This  correction  was  applied  during  MAPTIP  to  adjust  for  detector  response 
drift  (Chan,  1993). 

3.  Polarizers 

The  polarizers  utilized  during  MAPTIP  were  Cambridge  Physical  Sciences 
Model  IGP-228-71  with  a  diameter  of  71  mm.  The  polarizing  material  is  a  0.12  pm 
ainminnm  grid  of  0.25  pm  pitch  on  a  substrate  of  barium  fluoride.  They  have  85% 
transmissivity  and  98%  polarization  at  3.9  pm.  The  polarizers  were  placed  externally 
on  the  scanner  lens  systems  and  rotated  manually  to  achieve  vertical  and  horizontal 
polarizations.(Chan,  1993)  The  filters  placed  externally  reduce  the  system  sensitivity 
due  to  the  added  attenuation  effects  of  the  filters.  Any  change  in  system  setup,  such 
as  including  filters  requires  recalibration  on  the  AGA  780  (AGA  OP  MANUAL,  1980). 
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4.  AGEMA  CATS  Program 

The  AGEMA  CATS  software  allows  for  processing  images  immediately  or  for 
storing  images  to  be  analyzed  in  the  future.  The  program  is  used  by  two  PC’s,  each 
connected  to  one  of  the  scanner  channels  and  obtains  every  other  field  produced  by  the 
AGA  780.  Thus  CATS  generates  a  140x140  array  for  each  image  recorded. 

C.  IMAGE  PROCESSING  AND  ANALYSIS 

Although  processing  capability  is  available  with  AGEMA  CATS,  this  thesis 
utilizes  the  Interactive  Data  Language  for  all  image  processing  and  analysis. 

1.  Interactive  Data  Language  (IDL) 

The  Interactive  Data  Language  (IDL)  program  was  used  to  conduct  the  analysis 
in  this  thesis.  This  software  is  an  array-based  system  with  extensive  capabilities  in 
image  processing.  Numerous  analysis  methods  and  graphical  representations  are 
available  for  both  programming  and  interactive  manipulation.  All  images  and  most 
plots  displayed  in  this  thesis  were  produced  with  IDL  and  converted  to  TIFF  files.  The 
files  were  then  printed  with  Pagemaker4.  Three  IDL  programs  were  produced  in  the 
coiorse  of  this  thesis  and  are  presented  in  Appendix  A.  Also,  numerous  routines  were 
generated  to  help  with  the  analysis  of  images  and  subsequently  to  assist  in  the 
formulation  of  the  three  analysis  techniques  contained  in  the  programs  of  Appendix  A. 

2.  Image  Format 

As  stated  earlier,  all  images  were  recorded  by  AGEMA  CATS  for  processing  by 
IDL.  Each  image  is  stored  as  a  140x140  byte  file.  The  first  846  bytes  of  the  file 
contain  various  image  information  tags,  such  as  filename,  date,  time,  offset  correction 
applied,  etc.  An  example  of  this  "header"  information  is  shown  below. 


Thu  Nov  17  17:42:44  1994  Date  of  ANALYSIS 
e:\holland\octl9\lw 


filename 

level 

range  offset  over 

imder 

date 

time 

OCT19100.IMG 

48.80 

2.00  -32.00 

19600 

0 

10/19/93 

13:  4:54 

OCT19101.IMG 

44.20 

2.00  -32.00 

1226 

13379 

10/19/93 

13:13:11 

OCT19102.IMG 

44.20 

2.00  -32.00 

1175 

16313 

10/19/93 

13:13:19 

OCT19103.IMG 

44.40 

5.00  -32.00 

226 

5716 

10/19/93 

13:13:29 
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The  image  file  is  read  into  a  140x140  array  in  IDL.  The  array  values  are  the  scanner 
generated  thermal  values  scaled  from  0  to  255.  Zero  represents  the  coldest  value  in 
the  image  and  255  represents  the  hottest  value  viewed  (AGA  Manual,  1980). 

3.  Narcissus  Effect 

All  the  MAPTIP  images  show  a  dark  spot  in  the  center  of  the  image  similar  to 
that  shown  in  Figure  4.1.  This  is  known  as  a  Narcissus  spot.  The  Narcissus  spot  is 
created  by  retro-reflection  of  radiation  from  the  detector  plane  by  the  external  filters. 
Radiation  fi-om  the  cold  detector  mounting  at  approximately  77K  reflects  back  toward 
the  detector.  The  geometry  of  the  internal  optics  are  such  that  the  cold  detector 
radiation  focuses  on  the  detector  as  a  small  circular  region  near  the  center  of  the 
image.  The  effect  of  the  Narcissus  spot  in  addition  to  the  scene  obstruction 
encountered  is  to  reduce  the  effective  d5mamic  range.  The  extremely  low  values  of 
irradiance  from  the  Narcissus  relative  to  the  average  scene  value  compresses  the  range 
of  scene  values  in  the  upper  end  of  the  dynamic  range,  resulting  in  a  loss  of  scene 
resolution.(Lloyd,  1975) 

The  Narcissus  in  the  system  used  occurs  because  the  filters  are  moimted 
normal  to  the  radiation  from  the  detector  plane  collimated  by  the  lens  system.  The 
reflected  radiation  is  then  focused  back  on  the  detector  plane.  Without  the  filters 
installed,  the  Narcissus  is  not  present. 

Images  were  taken  with  regard  to  the  Narcissus  spot.  Target  images  were 
placed  away  from  the  central  Narcissus  spot  during  data  collection  to  reduce  its  affect. 
Although  the  Narcissus  could  in  principle  be  removed  by  image  processing  (Chan, 

1993),  it  would  be  difficult  if  not  impossible  to  replace  the  lost  information  accurately 
with  other  than  an  average  or  a  randomly  modulated  average  background  radiance 
value.  In  many  images,  a  substantial  depression  in  background  is  realized  that  is 
difficult  to  ascertain  from  a  strict  viewing  of  the  image  alone.  This  affect  can  be  seen 
in  Figure  4.11  as  not  only  the  hole  from  the  Narcissus,  but  also  the  depressed  trough 
spanning  the  length  of  the  image. 

4.  Ship  Location  Techniques 

Calculating  ship-to-background  contrast  requires  a  distinction  between  ship 
data  and  background  data.  The  ship  or  other  target  of  interest  must  be  located  within 
the  image  array  in  order  to  process  the  target  data  separately  from  the  backgroimd 
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data.  This  requires  a  contrast  between  the  ship  and  background  as  given  by  Equation 
3.41.  Three  primary  methods  for  processing  the  image  array  data  are  presented  below. 
All  three  utilize  a  thresholding  method  for  locating  the  ship.  Essentially,  thresholding 
means  acceptance  of  only  pixels  values  of  the  array  representing  brightness  that  are 
above  a  defined  threshold  and  setting  to  zero  all  values  below  this  threshold.  Ideally, 
all  ship  pixel  values  are  required  to  be  above  the  threshold  and  all  background  pixels 
values  below.  Ship  pixel  values  below  the  background  level  may  be  lost.  Technically, 
to  remove  all  the  background,  the  threshold  must  be  set  to  the  maximum  backgroimd 
value.  Loss  of  ship  radiance  values  below  this  value  is  possible. 

Care  must  be  taken  to  compare  appropriate  backgroimd  and  ship  values. 

Figures  4.1a  and  4.2a  show  tjpical  high  signal-to-noise  ratio  images  for  vertical  and 
horizontal  polarizations,  respectively.  Figures  4.1b  and  4.2b  depict  corresponding 
radiance  profiles.  The  dark  curves  represent  the  average  ship  radiance  profile.  This  is 
calculated  by  summing  the  ship  radiance  along  each  row  of  the  image  array  and 
dividing  by  the  number  of  ship  pixels  in  that  row.  The  gray  curves  depict  the 
background  radiance  profile.  This  is  calculated  below  the  horizon  by  averaging  all 
pixels  that  are  not  ship  pixels.  Above  the  horizon  the  average  is  calculated  using  only 
pixels  along  the  edge  of  the  image  to  avoid  the  Narcissus.  Figures  4.3a  and  4.4a  show 
similar  images  for  low  signal-to-noise  ratio  images.  Representative  background  profiles 
of  the  image  in  Figure  4.1a  are  depicted  in  Figure  4.5a  and  4.6a.  The  upper  curve  in 
Figure  4.5a  shows  the  sea  radiance  profile  near  the  bottom  of  the  image.  The  three 
lower  curves  show  the  sky  radiance  profiles.  The  topmost  of  the  three  curves 
represents  the  lowest  elevation.  The  lower  curves  represent  sequentially  increasing 
elevation.  Thus,  the  sky  radiance  decreases  with  elevation.  Figure  4.6a  shows  a 
vertical  slice  through  the  image  background.  Analysis  of  this  figure  and  Figure  4.5a 
reveals  that  the  background  varies  significantly  in  the  vertical  and  less  in  the 
horizontal  ignoring  effects  from  the  Narcissus  spot.  The  background  is  approximately 
constant  in  value  below  the  horizon,  but  rises  sharply  above  the  horizon  and  then 
tapers  off  at  higher  elevations.  Therefore,  to  accurately  define  the  contrast  between 
the  ship  and  background,  the  ship  values  below  the  horizon  should  be  compared  with 
background  values  below  the  horizon.  Likewise,  ship  values  above  the  horizon  should 
be  compared  with  background  values  above  the  horizon.  In  addition,  the  pixels  used  to 
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calculate  the  background  average  values  should  be  as  close  to  the  ship  as  possible  to 
represent  accurately  the  contrast  between  the  ship  and  the  background  in  the 
immediate  vicinity. 

Thresholding  works  well  for  targets  with  a  high  signal-to-noise  ratio.  However, 
choosing  an  appropriate  threshold  value  is  not  clearly  defined  and  somewhat 
subjective.  Figure  4.7a  shows  the  image  from  Figure  4.1a  and  4.2a  above  using  a 
threshold  value  equal  to  the  background  average  value  times  1.2.  As  the  threshold 
value  is  increased  to  (1.3)  x  (backgroimd  average)  in  Figure  4.7b  and  (1.5)  x 
(backgroimd  value)  in  Figure  4.7c  more  of  the  ship  is  deleted  along  the  edges.  From 
Figure  4.7a  it  is  clear  that  choosing  a  threshold  value  is  not  difficult  for  high  signal-to- 
noise  ratio  targets.  Figure  4.8  shows  the  low  signal-to-noise  ratio  image  of  Figures 
4.3a  and  4.4a  with  background  multiplying  factors  of  1.01,  1.03,  and  1.05,  respectively. 
In  this  series  of  images  some  background  exists  in  the  remaining  image  along  with  the 
ship.  In  this  particular  case,  the  backgroimd  can  be  separated  from  calculation  of  ship 
radiance  by  confining  the  calculations  to  a  region  near  the  ship.  However,  in  many 
images,  some  background  remains  too  close  to  the  ship  to  avoid  and  also  introduces  a 
varying  source  of  error.  As  the  background  factor  is  increased  sequentially,  thus 
increasing  the  threshold  values,  more  of  the  background  from  the  previous  image  is 
removed.  However,  increasing  amounts  of  the  ship  image  are  removed,  also. 
Appropriate  thresholding  is  a  balance  between  removing  a  sufficient  quantity  of 
background  while  retaining  a  representative  amount  of  the  target  ship.  This  balance 
becomes  increasingly  difficult  with  lower  signal-to-noise  ratio  targets.  Often,  a  large 
part  of  the  ship  is  not  visible  over  the  surrounding  background  and  cannot  be  located 
in  the  image.  These  concepts  are  further  described  and  applied  in  the  following  image 
processing  programs. 
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Figure  4.1  Vertical  Polarization,  for  high  signal-to-noise  ratio  image  (a)  IDL  genera,tecl  image  taken  at  1313 
Oct  19, 1993,  (b)  Average  ship  and  average  background  apparent  radiance  profiles.  Ship  radiance  is 
denoted  by  the  black  curve  and  background  radiance  by  the  gray  curve. 


Figure  4.2  Horizontal  Poiarization  for  high  signal-to-noise  ratio  image  (a)  IDL  generated  image  taken  al 
1313  Oct  19, 1993,  (b)  Average  ship  and  average  backaroond  apparent  radiance  profi,!es.  Ship  radiance 
denoted  bv  the  black  curve  and  background  radiance  by  the  gray  curve. 


Figure  4.3  Vertical  Polaiization  for  low  signal-to-noise  ratio  image  (a)  IDL  generated  image  taken  at  0915 
Oct  28, 1993,  (b)  Average  ship  and  average  background  apparent  radiance  profiles.  Ship  radiance  is 


denoted  by  the  black  curve  and  background  by  tlie  gray  curve. 


Figure  4.4  Hoiizontal  Polarization  for  low  signal-to-noise  ratio  image  (a) 
Oct  28, 1993,  (b)  Average  ship  and  average  background  apparent  radiance 
denoted  by  the  black  curve  and  background  radiance  by  the  gray  cur\'e. 


noise  ratio  image  (a)  Upper  curve  shows  sea 
background,  lower  three  curves  show  sky  back¬ 
ground,  with  highest  of  the  three  representing  the 
lowest  elevation(b)  Profile  through  ship.  Wings 
represent  sea  background. 
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Figure  4.6  Vertical  profiles  for  high  signal-to- 
noise  ratio  image  with  horizon  at  position  100  (a) 
Background,  (b)  Profile  through  ship.  Note  cold 
spot  from  Narcissus  where  radiance  drops  to  zero. 
Ship  located  at  large  spike  centered  on  position 
100. 
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Figure  4.7  Thresholding  high  S/N  ratio  image. 
Left  image  is  vertical  polarization  and. right  image 
is  horizontal  polarization  (a)  Background  factor  = 
1 .2  (b)  Background  factor  =  1 ,3  (c)  Background 
factor  =1.5. 


5.  Image  Processing  Programs 

This  section  presents  three  processing  programs  developed  to  calculate  the 
contrast  improvement  factor  associated  with  the  visual  improvement  in  contrast 
apparent  in  the  MAPTIP  images.  Although  contrast  improvement  in  horizontal 
polarization  over  vertical  polarization  is  visually  apparent  in  the  image  pairs,  the 
pursuit  of  these  programs  was  to  quantify  the  contrast  improvement  of  horizontal 
and/or  vertical  polarization  over  impolarized  images.  These  methods  are  limited  to  the 
apparent  radiance  received  and  recorded  as  image  files.  The  calculations  do  not 
consider  visual  modulation  transfer  functions  associated  with  the  detection  or 
identification  of  an  image  as  a  ship,  helicopter,  airplane,  etc.  The  computations 
accoimt  for  the  relative  differences  in  apparent  radiance  between  target  images  and 
backgroimd. 

Crucial  to  these  calculations  is  the  concept  of  generating  an  impolarized  image 
from  two  images  of  opposite  polarization.  Stokes  showed  that  light  intensities  of 
perpendicular  polarizations  can  be  added  to  obtain  an  unpolarized  image  or 
that  the  intensities  of  a  polarized  and  unpolarized  image  can  be  added  to  obtain  a 
partially  polarized  image  (Collett,  1993).  The  AGA  780  response  is  proportional  to  the 
intensity  of  radiation  received  (AGA  Operating  Manual,  1980).  Therefore,  the  addition 
of  a  vertically  polarized  image  with  the  same  image  in  horizontal  polarization  produces 
an  unpolarized  representation  of  the  image.  However,  the  CATS  program  stores  the 
images  scaled  relative  to  0-255  with  zero  being  the  lowest  radiance  received  and  255 
the  highest  radiance  received.  Also,  IDL  displays  images  as  bsde  array  values  on  the 
same  scale.  Therefore,  to  produce  the  unpolarized  image  with  IDL,  the  two  polarized 
images  are  averaged  as 


<N>+<N>,, 
<N>  = - ^ 

BP  O 


4.2 


Image  registration  is  another  important  element  in  the  process.  Since  the 
images  are  not  taken  simultaneously,  the  ship  has  moved  between  images.  Vertically 
this  is  usually  not  more  than  a  couple  of  scan  rows,  but  even  that  may  be  significant. 
However,  depending  upon  ship  speed  and  range,  the  movement  horizontally  may  be 
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substantial.  This  prevents  the  straight  addition  of  the  individual  arrays.  The  ships 
within  the  arrays  must  be  registered  on  to  each  other  in  order  to  average  their  values. 

CL  Contiguous  Area  Method 

The  contiguous  area  method  was  produced  by  W.J.  Lentz  at  the  Naval 
Postgraduate  School.  This  program  automatically  finds  the  hottest  spot  in  the  image 
and  works  through  a  contiguous  set  of  pixels  to  the  edges  of  the  ship  where  the 
temperature  drops.  The  spot  may  be  chosen  if  the  hottest  spot  is  backgroimd  clutter. 
Pixel  values  above  a  set  backgroimd  threshold  value  are  chosen  and  all  others  set  to 
zero.  The  background  threshold  value  may  be  varied.  Each  horizontal  line  in  the 
image  is  ascribed  a  background  temperature  based  on  a  smoothed  average  along  the 
edges  of  the  image.  The  values  are  used  to  find  variations  in  temperature  along  the 
horizontal  lines  that  could  represent  the  ship.  A  small  temperature  difference  might 
miss  much  of  the  ship.  For  this  reason  the  search  is  iterative  over  the  entire  image 
using  an  edge  threshold.  Once  the  ship  is  detected,  another  threshold  is  taken  in  a  box 
around  the  ship.  This  corrects  for  effects  such  as  the  Narcissus  spot  superimposed  on 
the  image.  The  output  is  a  file  that  gives  the  contrast  of  the  ship  to  background  in  a 
local  box  that  extends  to  the  limits  of  the  ship.(Lentz,  1994) 

The  thresholding  idea  was  first  used  to  discern  the  ship  from  the 
background  in  a  joint  experiment  with  NOSC.  In  that  experiment,  the  AGA  780  was 
used  without  polarizing  filters,  and  it  was  relatively  easy  to  find  the  ship.  In  the 
MAPTIP  experiment,  the  signal-to-noise  ratio  was  much  worse  both  because  of  the 
polarizers  that  reduced  the  signal  but  also  because  the  targets  were  generedly  much 
further  away.  The  earher  program  either  could  not  find  the  ship  or  would  include  too 
much  background  with  the  ship.  In  addition,  the  sink-hole  effect  in  the  center  of  the 
AGA  pictin-es  caused  by  the  Narcissus  effect  produced  additional  distortions  that 
caused  problems.  The  essential  idea  to  correct  these  problems  is  the  contiguous  point 
finding  program  in  IDL  called  search2D  used  in  the  subroutine  shipfind.pro.  It  starts 
from  the  designated  hot  spot  and  eliminates  noncontiguous  points  from  the  picture. 
Thus  what  remains  is  a  continuous  matrix  of  pixels  values  that  are  above  threshold. 

All  pixels  not  connected  to  the  hot  spot  are  set  to  zero.  Then  a  box  is  defined  at  the  x 
and  y  limits  of  the  ship.  Then  the  threshold  which  had  been  at  the  edge  of  the  picture 
is  redefined  for  the  vertical  edges  of  the  box  in  which  the  ship  is  contained.  The 
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process  is  repeated  producing  a  cleaner  ship  with  a  lower  allowable  temperature 
difference  for  backgroimd.(Lentz,  1994) 

Since  ship  pictures  are  taken  in  pairs,  the  program  takes  the  next  image 
and  compares  it  to  the  previous  image  and  combines  the  ship  images  to  produce  an 
impolarized  ship  image.  The  ship  may  have  moved  slightly  from  picture  to  picture,  so 
the  program  adjusts  for  the  motion  of  the  box.  One  polarization  or  the  other  will  have 
an  improved  contrast  over  the  backgrotmd  resvdting  in  a  larger  ship  image  with  more 
pixels.  In  some  cases,  for  example,  the  ship  will  consist  only  of  a  smokestack  with  one 
polarization  and  in  the  other  the  entire  ship  will  be  clearly  defined.  In  order  to 
properly  compare  the  contrasts,  the  ship  with  the  larger  number  of  pixels  is  used  to 
redefine  the  ship  outline  in  the  poorer  picture.  The  box  and  ship  sizes  are  the  same  for 
both  pictures  giving  a  better  calculation  of  the  polarization  improvement.  The  ship 
contrast  is  formed  between  all  the  points  on  the  ship  and  all  the  points  in  the  box 
containing  the  ship  that  are  not  the  ship.  This  is  done  for  both  polarizations  and  for 
the  sum  of  both  polarizations.  (Lentz,  1994) 

The  significant  advantage  of  this  program  is  its  autonomy.  Many  image 
pairs  can  be  processed  in  a  minimum  amount  of  time.  However,  this  program  does  not 
compare  the  entire  ship  as  observed.  Since  only  contiguous  areas  are  used,  only  the 
contiguous  area  containing  the  hot  spot  of  each  ship  image  is  used  to  calculate 
contrast.  h.  Outline  Method 

The  "outline  method"  is  an  attempt  to  overcome  the  difficulty  of  finding 
a  low  signal  to  noise  ratio  ship  image  and  is  provided  in  Appendix  A  as  the  program 
outlinlS.pro.  The  high  signal-to-noise  ratio  ship  image  of  Figure  4.1  was  used  to 
generate  a  block  outline  of  the  beam  aspect  of  the  ship.  The  outline  was  referenced  to 
an  IDL  generated  box  function.  This  allowed  for  processing  ship  targets  at  various 
ranges  since  at  longer  ranges  the  ship  takes  up  less  area  in  the  field  of  view.  When 
the  box  is  placed  around  the  ship  touching  only  the  stack,  bow,  stem,  and  waterline, 
the  outline  is  scaled  to  fit  the  ship.  If  only  part  of  the  ship  is  visible,  then  judgment  is 
required  to  align  the  box  properly. 
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Once  the  box  is  placed  around  the  ship,  average  background  values  are 
calculated  row  by  row  using  five  pixels  to  either  side  of  the  box.  The  outline  is  drawn 
referenced  to  the  box  as  shown  in  Figure  4.9  All  processing  was  performed  on  raw 
data  image  files.  Figure  9  is  an  expzmded  file  of  dimensions  280  x  280.  This 
expansion  process  consists  of  doubling  the  size  of  each  side  of  the  image  array.  The 
new  pixel  positions  are  filled  by  linearly  interpolating  between  the  original  pixels. 

Note  the  "steps"  created  during  the  expansion  along  the  bow  and  stem  outlines.  All 
pixels  outside  of  the  box  are  then  set  to  zero.  The  average  background  value  of  each 
row  multiplied  by  an  input  factor  to  account  for  pixels  above  the  average  is  used  to  set 
to  zero  any  remaining  pixels  within  the  outline  below  the  threshold.  Originally,  the 
backgroxmd  average  values  were  calculated  using  the  five  pixels  on  either  side  of  the 
outhne.  Ship  edge  effects  and  even  a  slight  error  in  placement  of  the  DDL  box  lead  to 
significant  errors  in  backgroimd  average  values  when  part  of  the  ship  registers  as  one 
or  more  of  the  10  backgroimd  pixels,  i.e.  part  of  the  ship  falls  outside  of  the  identified 
boundary.  Therefore,  five  pixels  on  either  side  of  the  box  were  used  to  generate 
average  background  values  for  each  row. 

At  this  point  all  that  remains  in  the  array  is  the  visible  portion  of  the 
ship.  Then  the  program  steps  through  each  row  calculating  the  average  ship  value  of 
that  row  using  the  remaining  pixels  in  that  row.  This  produces  a  profile  of  the  average 
ship  values  which  can  be  displayed  graphically  as  in  Figure  4.1b  and  4.2b.  A  second 
image  and  profiles  are  produced  for  the  perpendicular  polarization.  The  profiles  of 
contrast  and  contrast  improvement  may  also  be  represented  similarly.  Then  using 
these  values  average  contrast  and  average  contrast  improvement  for  above  and  below 
the  horizon  may  be  determined. 

However,  coupling  of  the  ship  areas  is  required.  Calculating  average 
ship  values  of  the  individual  polarizations  without  coupling  the  areas  is  comparable  to 
comparing  images  of  two  different  targets.  However,  this  method  erroneously 
accounted  for  this  coupling.  First,  the  "outline  method"  assumes  that  the  ship  images 
are  registered  exactly  due  to  the  outline  placement.  Therefore,  each  row  of  one  image 
is  compared  to  its  respective  row  in  the  orthogonal  polarization.  A  140  x  140  image 
window  on  the  computer  monitor  is  38mm  x  38mm.  This  requires  accurate  placement 
of  the  box  aroimd  the  ship  to  within  0.27mm.  The  assumption  of  proper  registration  is 
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invalid.  Second,  coupling  was  achieved  by  comparing  the  number  of  ship  pixels  in  a 
given  row  to  the  number  of  ship  pixels  in  the  same  row  of  orthogonal  polarization.  The 
row  with  the  most  pixels  and  thus  the  larger  area  must  represent  the  minimum  actual 
ship  area  along  that  elevation.  Then  the  sum  of  pixels  along  both  rows  was  divided  by 
the  larger  number  of  pixels  to  obtain  the  effective  average  ship  radiance  at  that  row. 
This  erroneously  reduces  the  apparent  average  radiance  of  the  row  with  the  lower 
number  of  pixels. 

Several  problems  exist  with  the  "outline  method"  in  addition  to  the 
coupling  described  above.  It  is  difficult  to  place  the  box  properly  aroimd  low  signal-to- 
noise  ratio  ship  images.  Improper  placement  and  sizing  results  in  the  outline  pattern 
not  including  all  of  the  ship  on  one  side  and  including  excess  background  on  the 
opposite  side.  The  major  drawback  with  the  "outline  method"  is  that  it  is  restricted  to 
beam  aspect  only.  The  image  array  may  be  inverted  without  loss  of  information  in 
order  to  utilize  the  outline  on  either  the  port  or  starboard  aspect.  However,  more 
detailed  programming  is  required  to  allow  for  other  than  beam  aspects  although  this 
method  is  probably  too  cmnbersome  to  provide  accmate  outlining  in  other  than  beam 
aspects. 

c.  Box  Method 

The  "box  method"  is  the  most  simple  of  the  three  methods  presented.  In 
fact  it  was  the  first  one  attempted,  but  later  dropped  to  permit  experimentation  with 
the  previous  methods.  This  method  simply  places  a  box  aroimd  the  ship  and 
determines  the  average  background  value  about  the  ship  row  by  row  using  the  five 
pixels  to  either  side  of  the  box.  All  pixels  outside  of  the  box  are  set  to  zero.  Then  a 
threshold  is  obtained  using  the  average  background  value  and  an  additional  amount  to 
account  for  backgroxmd  values  above  the  average.  As  before  the  values  below  this 
threshold  within  the  box  are  set  to  zero.  However,  this  method  contains  a  significant 
change  in  the  fundamental  interpretation  of  the  calculation  of  average  ship  value. 

The  "outline  method"  calculated  the  number  of  pixels  per  row  that  were 
ship  pixels  for  both  vertical  and  horizontal  polarizations.  The  larger  number  of  ship 
pixels  along  with  the  smn  of  ship  values  per  row  were  used  to  calculate  average  ship 
value  for  vertical  and  horizontal  polarizations.  The  two  polarizations  were  then 
averaged  to  produce  an  unpolarized  image.  This  method  was  flawed  because  it 
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improperly  couples  together  the  two  images  with  respect  to  ship  area.  The  image  with 
the  smaller  number  of  pixels  is  considered  to  have  an  area  greater  than  the  actual 
area.  Thus  pixel  values  of  zero  are  included  in  calculating  its  average  radiance.  Since 
the  ship  pixel  values  must  be  at  least  that  of  the  background,  assigmng  them  a  value 
of  zero  is  incorrect.  This  error  reduces  the  calculated  average  radiance.  The  "box 
method"  corrects  for  this  flaw. 

6.  Box  Method  Generation  of  Unpolarized  Image 

In  the  "box  method"  the  respective  ship  areas  in  the  vertical  and  horizontal 
polarizations  are  coupled.  The  ship  areas  are  represented  by  the  number  of  ship  pixels 
in  the  image.  Ship  area  is  subdivided  further  with  respect  to  ship  area  above  and 
below  the  horizon.  The  image  with  the  larger  area  or  greater  number  of  pixels  above 
or  below  the  horizon  defines  the  ship  area  for  that  region.  In  the  other  image  of 
smaller  area,  the  "missing"  pixels  are  given  a  value  equal  to  the  background  value. 

The  sum  of  all  pixels  in  their  respective  regions  is  divided  by  the  total  number  of  ship 
pixels  in  that  region.  Now  both  ship  images  contain  the  same  area  of  ship  surface. 

This  couples  the  ship  areas  by  referencing  the  average  ship  values  to  the  same  ship 
area.  The  motivation  for  this  is  that  while  part  of  the  ship  in  one  image  cannot  be 
seen  over  the  background  it  actually  is  being  seen  with  the  value  of  the  background. 
Now  the  average  ship  value  of  each  polarization  is  being  referenced  to  the  same  size  or 
area  of  ship  although  this  area  is  not  necessarily  the  true  ship  projected  area.  Even 
the  image  with  the  larger  area  may  not  represent  the  total  projected  area  if  part  of  its 
radiance  is  the  same  as  the  background  average.  Without  this  coupling  provided,  the 
smaller  ship  image  average  value  is  calculated  erroneously  lower  than  the  actual  value 
leading  to  erroneously  low  contrast  and  contrast  improvement  values.  In  certain 
instances,  the  "outline  method"  produced  negative  contrast  values  which  was 
impossible  since  the  detectable  ship  is  of  higher  radiance  than  the  background.  The 
negative  contrast  values  indicated  the  existence  of  and  led  to  the  identification  of  the 
area  coupling  flaw. 
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Some  image  information  is  lost  with  the  "box  method."  Since  the  average  ship 
values  above  and  below  the  horizon  are  calculated  as  the  sum  of  all  ship  radiance 
values  divided  by  the  total  number  of  pixels,  the  row  by  row  profile  information  is  lost. 
If  the  two  ship  images  could  be  accurately  registered,  then  a  row  by  row  comparison 
would  be  possible.  However,  since  the  relative  radiance  values  of  the  ship  images  with 
respect  to  polarization  do  not  differ  beyond  what  could  be  attributed  to  noise,  the 
profile  information  from  these  images  would  provide  no  new  meamngful  information. 

7.  Ship  Signature  Polarization 

In  an  attempt  to  analyze  the  polarization  differences  between  horizontal  and 
vertical  polarizations,  two  methods  were  used.  The  first  is  a  three  dimensional  surface 
plot  or  surface  representation.  The  second  is  a  cross-sectional  view  of  the  image  array 
or  profile  representation. 

a.  Siirfmx  R^TVsentcdwn 

A  three  dimensional  surface  plot  of  an  image  array  is  shown  in  Figure 
4.10.  The  top  image  is  of  the  vertical  polarization  of  a  ship  image  and  the  bottom  is  of 
horizontal  polarization.  Figure  4.11  shows  the  same  plot  rotated  about  the  vertical 
axis.  The  effect  of  the  Narcissus  can  be  seen  better  from  this  angle.  The  Narcissus 
can  be  seen  as  the  depression  or  hole  in  the  middle  of  the  plot.  Additional  damage  by 
Narcissus  depression  can  be  seen  emanating  away  from  the  primary  "hole"  in  the  form 
of  a  trough.  This  cold  spot  reduces  the  effective  dynamic  range  of  the  array  values  by 
forcing  the  target  radiance  values  to  be  compressed  at  the  upper  end  of  the  dynamic 

range. 

Comparison  of  the  two  surface  plots  shows  that  it  is  difficult  to  discern 
ship  radiance  value  differences.  What  differences  are  present  cannot  be  justified  as 
anything  other  than  noise.  Also,  it  is  difficult  to  determine  polarization  affects  on  the 
backgroimd.  Lastly,  this  data  provides  no  numerical  information,  although  it  does 
provide  some  visual  insight.  Some  images  such  as  the  Lynx  helicopter  in  Figure  4.12 
exhibit  flat  surfaces  in  these  plots  as  shown  in  Figure  4.13.  These  fiat  surfaces 
represent  pixels  values  that  exceed  255.  Their  actual  apparent  values  are 
undetermined. 
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b.  Profile  R^resentation 

An  alternate  representation  is  shown  in  Figure  4.5b  and  4.6b.  These 
show  horizontal  and  vertical  slices  through  the  ship.  All  the  slices  of  either  vertical  or 
horizontal  profiles  can  be  averaged  to  obtain  a  representative  profile  of  the  ship.  The 
Narcissus  spot  and  the  nature  of  the  background  lead  to  averaging  the  background 
pixel  value  of  each  row  instead  of  each  column  to  produce  a  representative  vertical 
profile.  Vertical  profiles  complement  the  process  of  comparing  ship  row  values  with 
their  respective  background  row  values.  Comparison  of  vertical  profiles  of  horizontal 
and  vertical  polarizations  again  shows  little  polarization  of  radiation  from  the  ship. 
However,  vertical  profiles  show  the  effect  polarization  filters  have  on  the  below-horizon 
backgrotmd  radiation  received.  This  is  demonstrated  in  Figures  4.14  and  4.15.  These 
graphs  show  superpositions  of  the  average  horizontal  and  vertical  backgrovmd  radiance 
values  and  ship  values.  The  impolarized  radiance  would  be  the  average  of  the  two  and 
positioned  midway  between  the  lines  shown.  The  background  radiation  received  below 
the  horizon  is  significantly  reduced  with  the  horizontal  polarization  filter.  A  small 
difference  in  radiation  received  exists  above  the  horizon.  Thus,  the  sea  backgrotmd 
radiation  is  predominantly  vertical  polarized.  No  discemable  difference  exists  along 
the  ships  in  Figures  4.14b  and  4.15b  and  thus  the  ship  radiance  is  essentially 
unpolarized. 

c.  Average  VcHue  CvdcuUdUm 

The  average  value  of  ship  radiance  for  a  given  polarization  may  be 
calculated  with  two  methods.  The  first  method  calculates  the  ship  average  values  for 
both  above  euid  below  the  horizon.  The  "box  method"  program  reduces  the  current 
image  array  to  only  ship  pixels  and  all  other  pixels  are  set  to  zero.  The  horizon  is 
marked  with  the  mouse  cimsor.  The  average  ship  value  above  the  horizon  is  calculated 
simply  by  summing  the  pixels  above  the  marked  horizon  and  dividing  by  the  number  of 
ship  pixels  above  the  horizon.  The  average  ship  value  below  the  horizon  is  calculated 
similarly. 
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Figure  4. 10  Surface  plot  of  image  in  Figure  4.  la.  Flatter  base  region  repre¬ 
sents  background  radiance  value.  Elevated  portions  show  ship  radiance 
values  (a)  Vertical  polarization  (b)  Horizontal  polarization. 


54 


Figure  4.1 1  Surface  plot  of  Figure  4.10  rotated  (a)  Vertical  polarization  (b) 
Horizontal  polarization. 
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5gure  4.12  IDL  generated  image  o: 


Figure  4.13  Surface  plot  of  helicopter  in  Figure  4.12.  Large  spike  shows  the  engine 
heating.  Smaller  peaks  are  from  local  hot  spots  on  the  helicopter  (a)  Vertical  polarization 
(b)  Horizontal  polarization. 
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Figure  4. 14  Characteristic  apparent  radiance 
profile  comparisons  for  high  S/N  ratio  image  (a) 
Background  comparison  (b)  Ship  comparison. 


Figure  4.15  Characteristic  apparent  radiance 
profile  comparisons  for  low  S/N  ratio  image  (a) 
Background  comparison  (b)  Ship  comparison. 


Alternatively,  the  average  value  of  ship  radiance  for  a  given  polarization 
is  calculated  for  each  row.  This  average  is  calculated  simply  by  sxmiming  each  row 
between  the  top  and  bottom  of  the  box  and  dividing  by  the  number  of  ship  pixels  in  a 
row.  These  averages  then  constitute  the  image  representative  vertical  ship  profile. 

8.  Background  Polarization 

Background  polarization  effects  were  analyzed  using  image  profiles.  Significant 
vertical  polarization  was  observed  below  the  horizon.  Less  vertical  polarization  was 
observed  above  the  horizon. 

a  Average  Value  Calculatum 

The  average  value  of  backgrotmd  radiation  for  a  given  polarization  is 
calculated  for  each  row.  The  procedure  varies  with  the  location  of  the  target  with 
respect  to  the  Narcissus  spot.  For  images  where  the  target  is  above  the  Narcissus 
spot,  the  backgroimd  above  the  ship  (i.e.  above  the  box)  is  calculated  simply  by 
summing  each  row  and  dividing  by  the  mnnber  of  pixels  in  a  row  (i.e.  140).  Below  the 
ship  the  first  and  last  20  pixels  in  each  row  are  summed  and  divided  by  40  to 
determine  the  average  background  value  below  the  ship.  Along  the  ship  the  five  pixels 
to  either  side  of  the  box  are  averaged  to  find  the  average  backgroimd  value  along  the 
ship.  These  averages  then  constitute  the  image  representative  vertical  background 
profile. 

6.  Characteristic  Profile 

Previously,  Figinre  4.5  and  4.6  showed  both  horizontal  and  vertical 
profiles  of  an  image.  The  horizontal  profiles  depict  a  backgroimd  that  is  essentially 
constant  at  a  given  elevation  disregarding  the  affect  of  the  Narcissus.  The  vertical 
profiles  depict  a  widely  varying  background  with  respect  to  elevation.  A  target  is 
discemable  from  the  background  in  the  immediate  vicinity  if  there  exists  an 
appreciable  contrast  between  the  ship  radiance  and  the  background  radiance.  The  key 
is  "immediate  vicinity."  Contrast  must  be  based  upon  appropriately  related  ship  and 
background  values.  It  is  not  appropriate  to  compare  ship  radiance  values  below  the 
horizon  with  background  above  the  horizon.  Therefore,  the  ship  radiance  values  along 
a  given  row  must  be  compared  with  background  radiance  values  along  that  same  row. 
Thus,  all  the  pixels  in  each  row  are  averaged  to  generate  a  characteristic  profile. 
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The  characteristic  profiles  of  Figures  4.14a  and  4.15a  show  a  significant 
reduction  in  background  radiance  below  the  horizon  while  no  significant  difference 
exists  in  the  radiation  from  the  target  ship.  It  is  this  reduction  in  background 
radiation  that  polarization  filters  can  take  advantage  of  and  increase  target  to 
background  contrast.  By  aligning  the  filters  to  block  the  predominantly  vertically 
polarized  background  radiation  the  contrast  will  be  enhanced. 

9.  Ship  to  Background  Contrast 

Once  the  ship  and  background  average  values  have  been  determined,  the 
contrast  values  may  be  calculated.  Contrast  is  defined  as  a  difference  in  radiance 
value  between  the  target  and  surrounding  backgroimd.  This  difference  may  be  positive 
or  negative  depending  on  whether  the  target  is  of  greater  or  lesser  radiance  than  the 
background.  Equation  4.3  below  provides  the  mathematical  definition  of  contrast. 
Ideally,  contrast  results  from  any  one  or  a  few  adjacent  pixels  being  of  greater  or  lesser 
value  tbfm  the  svurounding  backgroimd  pixels.  This  introduces  the  concept  of  clutt  -r. 
Clouds  may  reflect  radiation  into  the  detector  with  a  value  comparable  to  the  ship 
radiance.  The  ship  targets  contain  more  than  one  individual  pixel  as  is  necessary  to 
achieve  classification  or  identification.  So  there  exists  a  question  as  to  what  is 
detection  that  is  not  addressed  by  this  thesis.  As  stated  before,  this  thesis  does  not 
include  modulation  transfer  functions  for  determining  detection. 

10.  Contrast  Improvement 

The  quantities  of  interest  are  the  improvement  in  vertical  or  horizontal 
polarization  contrast  over  the  unpolarized  contrast.  Provided  there  is  a  significant 
improvement  in  one  of  the  polarizations,  then  such  a  polarizer  may  prove  useful  in  an 
actual  target  detection  or  tracking  system.  It  is  reasonable  that  a  variable  plane  of 
polarization  filter  might  be  included  in  system  design. 

a.  Defhntionof  Contrast  Improvement 

Target  to  backgrotmd  contrast  can  be  defined  as  the  normalized 

^tgt~^bkgd  4  3 

^bkgd 

radiance  difference  between  target  and  backgroimd  where  Njgt  and  N^kg^  axe  apparent 
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radiance  values,  respectively.  Therefore,  the  contrast  for  vertical,  horizontal,  and 
unpolarized  target  images  is 


(^vW 


4.4 


^  _^Wgt  ^^l)bkgd 


45 


Then  the  contrast  improvement  for  vertical  and  horizontal 

^KP- 

b.  Background  Dependency 
As  stated  before,  contrast  calculations  are  dependent  upon  the 
backgroimd  utilized  in  the  calculations.  If  the  backgroimd  average  is  determined  using 
the  entire  range  of  values  from  the  lowest  radiance  at  high  elevations  to  the  highest 
radiance  just  above  the  horizon,  then  this  average  may  be  significantly  higher  or  lower 
fVian  the  average  background  value  in  the  vicinity  of  the  target.  Care  must  be  taken  to 
associate  the  appropriate  background  with  the  target.  This  is  most  effectively  done  by 
comparing  the  below  horizon  background  with  the  ship  component  below  the  horizon 
and  the  above  horizon  background  with  the  ship  component  above  the  horizon.  A  more 
accurate  method  would  be  to  do  a  row  by  row  comparison.  This  is  easily  done  with  one 
image  of  vertical  or  horizontal  polarization.  However,  for  meaningful  analysis  the 
images  must  be  summed.  This  leads  to  the  problem  of  registration  covered  earher. 


polarizations  is  defined  as 
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a  Ship  Area  Reference  Dependency 

The  contrast  calculations  depend  on  the  area  coupling  discussed 
previously.  Originally,  row  by  row  analysis  was  performed.  This  entailed  summing 
the  ship  pixels  along  each  row.  The  same  size  box  was  used  in  each  polarization  and 
assumed  to  register  the  ship  image  properly.  The  same  row  of  each  image  was 
compared  by  using  the  larger  number  of  ship  pixels  in  each  row  from  the  respective 
vertical  and  horizontal  images.  This  often  leads  to  negative  contrast  values,  since  the 
calculated  N^gt  may  be  erroneously  low.  This  problem  was  resolved  by  first 
acknowledging  the  error  in  ship  registration.  Then  instead  of  a  row  by  row  analysis 
the  ship  radiance  was  calculated  as  an  average  either  above  or  below  the  horizon. 

Also,  the  difference  in  ship  pixels  was  compensated  for  by  assigning  the  background 
average  to  the  "missing"  pixels. 

D.  RESULTS  OF  ANALYSIS 

Contrast  improvement  calculations  were  generated  by  the  "contiguous  area" 
and  "box"  methods.  Visual  review  of  the  MAPTIP  image  data  reveals  a  greater  ship  to 
background  contrast  in  the  horizontally  polarized  images  than  in  the  vertically 
polarized  images.  Of  interest  is  the  degree  of  contrast  improvement  achieved  with 
polarizers  over  that  without  polarization  filters.  Unpolarized  image  data  is  not 
available  for  these  images  and  must  be  artificially  constructed.  Both  the  contiguous 
area  and  box  methods  generate  unpolarized  images  and  calculate  the  contrast 
improvement  obtained  with  vertical  and  horizontal  filters  as  described  above.  The 
results  of  these  calculations  are  provided  in  Appendix  B. 

Both  methods  show  cont  ast  improvement  in  the  horizontal  polarization. 
Sometimes  even  vertical  polarization  filtering  shows  a  contrast  improvement  over 
xmpolarized  mages.  Horizontal  polarization  consistently  provides  a  greater  contrast 
improvement  over  vertical  polarization  except  within  a  strong  sunglint  channel.. 

Tables  4.1  through  4.5  show  results  from  selected  images  throughout  different 
times  of  the  day.  Table  4.1  images  were  recorded  in  the  early  morning.  Table  4.2  in 
mid-moming,  Table  4.3  in  mid-day,  Table  4.4  in  mid-aftemoon,  and  Table  4.5  in  the 
evening.  Each  table  is  in  two  parts.  Part  (a)  lists  the  average  radiance  values  for 
target  and  background  both  above  and  below  the  horizon  as  calculated  from  the  "box 
method."  Part  (b)  lists  contrast  and  contrast  improvement  values.  The  last  three 
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columns  depict  polarized  target-to-background  contrast  (C),  unpolarized  target-to- 
background  contrast  (Cup),  and  contrast  improvement  (Cj^jp)  as  calculated  from  the 
"contiguous  area  method." 

These  values  of  contrast  improvement  require  qualification.  The  results  of  both 
methods  are  not  strictly  comparable.  Errors  are  inherent  in  both  methods  due  to 
misidentification  of  the  horizon  or  ignoring  horizon  effects.  The  contiguous  area 
method  calculates  a  ship  contrast  irrespective  of  the  horizon.  This  is  not  a  significant 
miscalculation  in  most  low  signal-to-noise  ratio  images.  In  these  cases,  the  horizon  is 
seldom  visible  and  the  apparent  radiance  just  above  and  below  the  horizon  is  similar. 
This  can  be  seen  from  the  example  low  signal-to-noise  ratio  image  of  Figure  4.4.  The 
results  of  calculations  for  this  image  are  tabulated  in  Table  4.2  as  file  OC28L025  and 
OC28L026.  The  "box  method"  contrast  and  contrast  improvement  below  the  horizon  is 
nearly  the  same  as  the  overall  contrast  and  contrast  improvement  calculated  from  the 
"contiguous  area  method."  However,  in  images  where  the  horizon  difference  is  large, 
this  error  increases  substantially.  Image  files  OC19L105  and  OC19L106  of  Table  4.3 
show  the  results  for  the  example  high  signal-to-noise  ratio  image  of  Figure  4.1.  Since 
in  such  images  the  sky  radiance  just  above  the  horizon  is  much  larger  than  below  the 
horizon,  then  the  overall  ship  to  backgroimd  contrast  is  falsely  reduced  because  the 
average  backgroimd  is  erroneously  increased  by  averaging  above  and  below 
background  values.  This  results  in  the  negative  contrast  values  of  several  image  pairs 
in  the  "contiguous  area"  method  data  as  seen  in  Table  4.2b  for  file  OC25L084.  The 
"box  method"  avoids  this  error.  However,  accurate  identification  of  the  horizon  is 
difficult.  This  also  leads  to  an  error  in  calculation  of  the  above  and  below  horizon 
background  values. 

Table  4.6  compares  the  "box  method"  contrast  improvement  (Cg)  for  horizontal 
polarization  below  the  horizon  and  "contiguous  area  method"  for  contrast  improvement 
(Cg)  for  horizontal  polarization  to  two  previous  efforts  to  calculate  contrast.  The 
contrast  in  was  calculated  without  correctly  coupling  the  ship  areas.  In  this  method 
the  ship  image  of  smaller  area  was  assumed  to  have  the  same  area  as  the  larger  ship 
image.  The  total  pixel  radiance  value  was  divided  by  this  larger  area  which  means  the 
missing  pixels  were  set  to  zero.  This  reduces  the  average  ship  radiance  and  contrast 
for  that  polarization.  This  reduced  ship  radiance  also  reduces  the  unpolarized  ship 
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radiance  and  contrast  which  is  calculated  from  both  polarized  radiance  values. 
Therefore,  opposite  polarization  contrast  improvement  will  be  increased  when  its 
contrast  is  divided  by  the  reduced  impolaiized  contrast.  If  the  vertically  polarized 
image  is  of  less  area,  then  the  horizontal  contrast  improvement  will  be  inflated  as 
shown  in  Table  4.6  as  compared  to  C2  which  sets  the  missing  pixels  values  of  the 
smaller  image  equal  to  the  background  radiance.  The  last  column  shows  ship-to- 
background  contrast  calculated  from  mean  radiance  data  shown  by  Chan  (Chan,  1993). 
This  mean  radiance  data  consists  of  the  average  ship  or  background  radiance  within  a 
box  placed  either  within  the  ship  target  or  a  box  containing  background  only. 

The  "contiguous  area"  method  is  semi-autonomous.  Once  a  threshold  value  is 
input,  the  only  operator  action  necessary  occms  when  the  program  experiences  a 
conflict  in  location  of  the  ship  hot  spot.  This  spot  is  identified  by  the  operator  and  the 
program  continues.  In  several  instances,  the  "contiguous  area"  method  includes  a 
large  and  significant  portion  of  background  in  its  calculation  of  ship  average  radiance 
value.  The  "box  method"  also  avoids  this  by  allowing  the  operator  to  set  the  size  of  the 
box  to  a  tight  fit  aroxmd  the  ship. 

Also,  the  "contiguous  area"  method  occasionally  does  not  include  all  of  the  ship 
in  calculating  average  radiance.  This  happens  when  portions  of  the  ship  in  both 
horizontal  and  vertical  polarizations  are  in  separate  blocks  of  pixels.  For,  example 
some  images  will  have  the  ship  stack  appear  separated  from  a  larger  area  of  ship  huU. 
The  program  will  only  use  one  of  these  areas  in  the  calculation.  For  strict  contrast, 
the  stack  area  which  is  of  a  greater  apparent  temperatm-e  should  be  chosen.  This 
leads  again  to  the  concept  of  detection.  If  the  stack  is  to  be  chosen  only  because  of 
higher  apparent  temperature,  then  in  all  images  only  the  highest  ship  pixel  value 
should  be  compared  to  the  surrounding  pixels.  For  the  stack,  the  surroimding  pixels 
are  other  ship  pixels.  The  concept  of  pixel  to  pixel  contrast  works  if  the  target  is  one 
pixel  in  area.  The  "contiguous  area"  method  compares  the  ship  "hot  chunk"  to 
backgromid  whereas  the  "box"  method  compares  the  entire  infrared  visible  ship  to 
background.  The  former  applies  more  to  pixel  detection  and  the  latter  applies  more  to 
target  detection. 

Apparent  in  the  "box"  method  data  is  the  greater  contrast  improvement  in 
below  the  horizon  compared  to  that  above  the  horizon  in  horizontal  polarization.  This 
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results  from  the  sea  background  emission  which  is  predominantly  vertically  polarized 
(Cooper,  Crittenden,  Milne,  Walker,  Moss,  and  Gregoris,  1992).  The  sky  background 
immediately  above  the  horizon  is  predominantly  unpolarized  where  there  exists  a 
significant  horizon  contrast  as  depicted  in  Figure  14a.  This  figure  shows  no 
measurable  change  in  apparent  radiance  in  either  polarization.  This  can  only  occur  if 
the  radiance  is  impolarized. 

Definitions  of  the  results  in  Tables  4.1  through  4.5  are  given  below.  The 
following  values  were  calculated  by  the  "box  method." 

(N)abv  =  apparent  radiance  above  the  horizon  for  given  polarization 

(N)bei  =  apparent  radiance  below  the  horizon  for  given  polarization 

(Nup)abv  =  apparent  radiance  above  the  horizon  for  unpolarized  image 

(Nup)bei  =  apparent  radiance  below  the  horizon  for  unpolarized  image 

Cabv  =  contrast  above  the  horizon  for  given  polarization 

Cbei  =  contrast  below  the  horizon  for  given  polarization 

(Cup)abv  =  contrast  above  the  horizon  for  unpolarized  image 

(Cup)bei  =  contrast  below  the  horizon  for  impolarized  image 

(Cbnp)abv  =  contrast  improvement  above  the  horizon  for  given  polarization 

(Cimplbei  =  contrast  improvement  below  the  horizon  for  given  polarization 

The  following  values  were  calculated  by  the  "contiguous  area  method." 

C  =  contrast  for  given  polarization  without  regard  for  horizon 
C„p  =  contrast  for  impolarized  image  without  regard  for  horizon 
Cjmp  =  contrast  improvement  without  regard  for  horizon 
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lable  4.1a  Kadiance  values  tor  selected  vertical/horizontal  image  pairs  and  unpolarized  image  with  respect  to  target  or 
background  source  and  horizon.  Images  taken  during  the  early  morning.  "V"  and  "H"  denote  vertical  and  horizontal  polarizati 
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Table  4.3a  Radiance  values  for  selected  vertical/horizontal  image  pairs  and  unpolarized  image  with  respect  to  target  or  ' 

background  source  and  horizon.  Images  taken  during  the  mid-morning.  "V"  and  "H"  denote  vertical  and  horizontal  polarizations, 
respectively. 
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Table  4.5a  Radiance  values  for  selected  vertical/horizontal  image  pairs  and  unpolarized  image  with  respect  to  target  or 
background  source  and  horizon.  Images  taken  during  the  mid-morning.  "V"  and  "H"  denote  vertical  and  horizontal  polarizations, 
respectively. 
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Table  4.5b.  The  first  six  columns  show  coi 
horizontal  polarization  as  calculated  by  "b 

FILE 

OC22L235V 

OC22L236H 

OC22L237V 

OC22L238H 

OC22L239V 

OC22L240H 

OC22L243V 

OC22L244H 

OC22L245V 

OC22L246H 

TIME 

20:58:50 

21:01:26 

21:02:07 

21:03:24 

21:03:43 
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calculated  by  "contiguous  area  method."  Images  taken  during  the  mid-morning.  "V"  and  "H"  denote  vertical  and  horizontal 
polarizations,  respectively. 


Image 

Pair 

File  Name 

Pol. 

Cl 

C2 

C3 

C4 

1 

OC19L087/88 

V/H 

7.7 

1.33 

1.00 

1.15 

2 

OC19L103/04 

V/H 

2.8 

1.74 

1.17 

1.14 

3 

OC19L105/06 

V/H 

1.33 

1.33 

1.16 

1.14 

4 

OC19L107/08 

V/H 

1.24 

1.23 

1.10 

1.13 

5 

OC19L109/10 

V/H 

1.22 

1.15 

1.12 

1.14 

6 

OC19L111/12 

V/H 

1.21 

1.13 

0.97 

1.11 

7 

OC19L122/23 

V/H 

1.24 

1.14 

1.07 

1.12 

8 

OC19L124/25 

V/H 

1.24 

1.14 

1.06 

1.10 

Table  4.6  Contrast  improvement  comparison.  calculated  without  setting  "missing" 


pixels  in  smaller  image  to  backgroimd  value,  Cg  calculated  "box  method",  C3  calculated 
using  "contiguous  area  method,"  C4  calculated  using  mean  radiance  data  obtained  by 
Chan.  Cl  and  Cg  represent  contrast  improvements  below  the  horizon  for  horizontally 
polarized  images.  C3  and  C4  are  contrast  improvements  for  horizontal  polarizations. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  below  horizon  contrast  improvement  values  confirm  that  a  significant 
improvement  in  detection  ability  is  possible  with  horizontal  filtering.  The  below 
horizon  values  are  significant  because  in  most  operational  environments  a  sea 
background  will  dominate  target  scenes.  For  an  aircraft  flsdng  as  low  as  500  feet,  the 
distance  to  the  horizon  is  25  nautical  miles  given  by 

D=1.14V(/i) 

where  D  =  distance  to  horizon  in  nautical  miles 
h  =  height  of  observer  in  feet 

(Tobin,  1974).  Then  the  distance  to  a  150  ft  high  surface  vessel  that  is  completely 
below  the  horizon  is  22.5  nautical  miles  as  shown  below. 


Figure  5.1  Observer  to  horizon  geometry. 


The  sxxrface  is  approximated  to  be  flat  in  the  diagram.  The  target  ship  height  is 
denoted  by  TS,  the  observer  height  is  denoted  by  OP,  and  the  horizon  by  H.  FLIR 
acquisition  ranges  are  often  displayed  as  a  function  of  target-to-backgroimd 
temperature  difference.  For  a  typical  400  ft  long  destroyer  from  a  beam  aspect, 
imclassified  acquisition  ranges  show  a  maximum  acquisition  range  of  less  than  20 
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nautical  miles  (Cooper  class  notes,  1994).  Thus  a  significant  operational  improvement 
is  possible  with  infi-ared  detection  and  imaging  systems  using  horizontal  polarizing 
filters  because  the  filters  would  reduce  sea  component  of  background  radiance. 

These  calculations  are  for  detection  only.  The  "box"  method  more  closely 
approximates  recognition  improvement  by  using  the  entire  ship  radiance  in  the 
average  radiance  calculation.  Further  consideration  needs  to  be  applied  to  the 
resolution  requirements  necessary  to  recognize  a  detected  target.  Target  detection  in 
the  presence  of  clutter  is  felt  by  some  to  require  some  resolution  capability  -  e.g.  FLIR 
resolution  requirements  are  established  at  1.5  line  pairs  for  detection  (Rosell  and 
Harvey,  1979).  Given  a  system  resolution  capability  requirement,  the  improvements  in 
detection  and  recognition  of  horizontal  filtering  over  no  filtering  should  be  explored. 

As  an  example,  some  image  pairs  show  a  sizeable  increase  in  total  ship  area  detected 
with  horizontal  polarization  (see  Figure  4.9).  A  system's  resolution  requirements 
would  be  less  with  the  increased  area  of  ship  radiance  received  in  excess  of  backgroimd 
radiance  received. 

It  is  difficult  to  develop  a  technique  for  evaluating  the  contrast  which  would  be 
observed  and  be  applicable  to  detection/recognition  criterion.  As  a  TninimnTn  the 
method  chosen  must  be  consistent  with  visual  results.  The  "box  method"  is  believed  to 
be  such  a  method.  Validation  of  that  conclusion  is  needed.  Contrast  improvement 
refers  only  to  apparent  contrast.  This  must  be  dependent  on  the  range  at  which  the 
target  is  viewed.  Range  dependence  of  the  contrast  should  be  included  and  validated. 

Image  registration  remains  a  problem.  The  inability  to  properly  sum  two 
images  precludes  correct  generation  of  an  unpolarized  image.  A  related  problem  is  the 
time  difference  between  image  pairs.  The  atmosphere  is  d3mamic  and  thus  the 
radiance  received  is  time  dependent.  The  atmospheric  changes  between  images  are 
unknown  and  may  or  may  not  be  significant. 

One  possible  solution  to  both  problems  is  in  progress.  This  involves  a  lens 
adapter  system  to  be  mounted  on  the  AGA  780.  This  system  serves  two  purposes. 

First,  the  design  is  such  that  the  Narcissus  spot  wdll  be  removed.  Second,  the  system 
splits  the  incoming  beam  of  radiation.  Each  beam  then  passes  through  either  a 
horizontal  or  vertical  polarizing  filter.  The  resulting  image  will  contain  both 
polarizations  each  140  x  70.  Since  these  images  were  taken  simultaneously,  t.bpn  they 
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will  be  registered  perfectly.  Then  they  may  be  summed  to  produce  an  unpolarized 
image.  Once  this  system  is  constructed  and  mounted  on  the  AGA  780  it  must  be 
calibrated  and  tested  to  verily  or  disprove  the  assertion  that  an  unpolarized  image  may 
be  produced  with  the  AGA  780  by  summing  a  vertical  and  horizontal  image. 
(Skowroneck,  1994) 

Environmental  data  acquired  from  the  MAPTIP  experiment  could  be  used  for 
further  analysis  of  the  MAPTIP  images.  Calculations  in  this  thesis  were  based  on 
apparent  radiance.  This  is  the  radiance  actually  received  at  the  AGA  780.  This  is  not 
the  radiance  emitted  from  the  target.  Contrast  and  contrast  enhancement  are  relative 
-  based  on  relative  radiances.  This  does  not  pose  a  problem  as  would  absolute 
radiances.  With  the  environmental  data  from  MAPTIP,  an  atmospheric  transmission 
model  such  as  modified  LOWTRANO  may  be  used  to  determine  the  actual  radiance  of 
the  ship  and  backgroimd.  This  information  is  necessary  for  proper  modelling  of 
electro-optic  system  requirements. 
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APPENDIX  A. 


A-  OUTLINE  METHOD  PROGRAM 


;outlinl4  constructs  outline  of  ship  referenced  to  box  placed  around  ship 
box  must  touch  top,  bow  tip,  stem  tip,  waterline 

!order=0 

print,'flirs  makes  movie  of  Airs' 
close,  1 

;close  in  case  of  crash 

intAi=0 

shpnmbr=l 

ct=0 

loadct,ct 

;window,0,xsize=280,ysize=280,xpos=0,3rpos=0 

on_ioerror,problem 

pathl='  ' 

Alel=' 

print, 'Ale  size  must  be  20446  with  header' 
read,'input  path  to  *.img  Ales',pathl 
if  pathl  ne  "  then  path=pathl 
print,'path=  ',path 

read,'input  type  of  Ales  eg  *.img  or  names  to  search  for',Alel 

if  Alel  eq  "  then  Alel='*.img 

Ales=AndAle(path+'  \  '+Alel,covmt=num) 

for  i=0,num-l  do  print,Ales(i),i 

loop:  print,n_elements(Ales)-l,'  last  number' 

read,'input  i  for  number  to  plot  or  -1  to  stop',i 

if  i  It  0  then  goto,  problem 

istart=i 

read,'mput  stop  number' ,istop 

for  i=istart,  istop  do  begin 

f=Ales(i) 

openr,l,f 

print,'opening  ',f,i 

head=bytarr(846) 

y=b3d;arr(140,70) 

yl=b3i:arr(140,70) 

37y=b5d;arr( 140, 140) 
readu,l,head 
readu,l,y 
readu,l,yl 
close,  1 

index=indgen(7  0) 

37y(*,index*2)=yl(*,index) 
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yy(*  ,index*2+ l)=y(* , index) 
ifl^maxCyy)  eq  255)  then  begin 
yl=yy(  where(yy  gt  254)) 
nover=n_elements(yl) 
endif  else  begin 
nover=0 
endelse 

jry=yy(*,  139-indgen(  140)) 

;save  ship  images 

.***^=**********************=1=*********** 

shipl=bytarr( 140, 140) 
ship2=bjd;arr(140, 140) 
if  shpnmbr  eq  1  then  begin 
shipl=yy 
shipll=yy 
shpmnbr=2 
titlel='vert' 

window, l,xpos=0,ypos=0,xsize=140,ysize= 140  ,title=titlel 
tvscl,shipl 
endif  else  begin 
ship2=yy 
title2='hor' 

window,2  ,xpos=0 ,3q)os= 160,xsize= 140,ysize=140  ,title=title2 
tvscl,ship2 
shpnmbr=l 
endelse 

endfor 

.***Sc************************************** 

;invert  outline  for  port  beam  aspect 

.*****:ic***************************^:*****H:** 

read,'invert  outline  for  port  beam  aspect?  yes=l  no=0', aspect 
if  aspect  eq  1  then  begin 
shipll=shipll(139-indgen(140),*) 
ship2=ship2(139-indgen(140),*) 
wset,l 
tvscl,shipll 
wset,2 
tvscl,ship2 
endif 

.***********jjjygj.|;  gjgjjj^J********** 

yl=l 

y=i 

;*********reset  variables********** 
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yy=congrid(yy , 140*2 ,140*2 ,Interp= 1) 

;!order=l 

loadct,0 

;wset,0 

print, 'Scaled  plot  of  data  in  O' 

;tvscl,yy 

range=double(head,770) 

level=double(head,762 ) 

dlevel=double(head,618) 

emiss=double(head,626) 

ct=double(head,650) 

r=double(head,570) 

b=double(head,5 78) 

f=double(head,586) 

anit=double(head,674) 

at=double(head,666) 

pl=emiss*ct 

p2=(  l-emiss)/emiss 

p3=1.0-ct 

ianib=r/(exp(b/amt)-f) 

iatm=r/(exp(b/at)-f) 

y=range*(float(yy)-128.)/256.0+level+dlevel 
y=y/pl-p2*iamb-p3/pl*iatm 
t=b/alog(abs((r/y+f)))-273.15 
;  get  rid  of  old  variables 
;delvar,y,yl,yy,liead 

shpave  l=fltarr(  140) 

shpavel(*)=0.0 

shpave2=fltarr(140) 

shpave2(*)=0.0 

shptotl=fltarr(140) 

shptotl(*)=0.0 

shptot2=fltarr(  140) 

shptot2(*)=0.0 

;indexl=fltarr(  140) 

;indexl(*)=0 
;index2=fltarr(  140) 

;index2(*)=0 
shppixl=fltarr(  140) 
shppixl(*)=0.0 
shppix2=fltarr(  140) 
shppix2(*)=0.0 
outline=fltarr  ( 140 , 140 ) 
outline(*,*)=0.0 
outline  l=fltarr( 140, 140) 
outlinel(*,*)=0.0 
outline2=fltarr( 140 , 140) 
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outline2(*,*)=0.0 


redo:read, 'enter  BF  value', dd 

BF=dd 

image=l 

nextimg: 

wset, image 

zz=boxfun(xO,yO,nx,ny,init=intfh/message) 

intfii=l 

box=tvrd(xO,yO,nx,ny) 

bkgdave=fltarr(  140) 

bkgdave(*)=0.0 

deltal=0 

delta2=l 

.*********************************************************** 

> 

;box  or  outline 

.*********************************************************** 


J 

;assign  image  array  to  outline  array  and  mark  horizon 

.*********************************************************** 

> 

if  image  eq  1  then  begin 
outline=shipll 
boat=shipll 

print, 'click  on  horizon  or  bottom  to  exit' 
replotl:  cursor,xv,yv,3 /device 
if  yv  It  5  then  goto,exitl 
xvl=xv 
yvl=yv 

print,'horizon=',yvl 
print,'chck  on  point  or  bottom  to  exit' 
goto, replotl 
exitl;print,'exiting' 
endif 

if  image  eq  2  then  begin 
outline=ship2 
boat=ship2 

print,'click  on  horizon  or  bottom  to  exit' 
replot2:  cursor,xv,yv,3/device 
if  yv  It  5  then  goto,exit2 
xv2=xv 
yv2=yv 

print, 'horizon=',yv2 
prmt,'click  on  point  or  bottom  to  exit' 
goto,replot2 
exit2  :print, 'exiting' 
endif 

.*********************************************************** 
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;calculate  background  average  values 

if  yO+ny/2.0  It  69  then  begin 
;box  below  narcissus 
;background  average  below  box 
for  n=0,y0-l  do  begin 
bkgdave(n)=total(outline(0:139,n))/140.0 
endfor 

;background  average  around  box 
for  n=yO,yO+ny-l  do  begin 

bkgdave(n)=(total(outline(x0-5:x0-l,n))+total(outline(x0+nx:x0+nx+4)))/10.0 

endfor 

;background  average  above  box 
for  n=y0+ny,139  do  begin 

bkgdave(n)=(total(outline(0;19,n))+total(outline(120:139,n)))/40.0 

endfor 

endif  else  begin 
;box  above  narcissus 
;background  average  below  box 
for  n=0,y0-l  do  begin 

bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 

endfor 

;background  average  around  box 
for  n=yO,yO+ny-l  do  begin 

bkgdave(n)=(total(outline(x0-5:x0-l,n))+total(outline(x0+nx:x0+nx+4)))/10.0 

endfor 

;background  average  above  box 
for  n=y0+ny,139  do  begin 
bkgdave(n)=(total(outline(0:139,n)))/140.0 
endfor 
endelse 

5 

;generate  outline  and  ship  average  values 

outline!*, 0:34)-1)=0.0 
;waterline 

templ=x0+(9.0/94.0)*nx 
temp2=x0+(89.0/94.0)*nx 
for  m=templ,temp2  do  begin 

if  outline(m,yO)  le  BF*bkgdave(yO)  then  outline(m,y0)=0.0 
endfor 

boat(templ:temp2,y0)=255 
outline(0:templ-l,y0)=0.0 
outline(temp2+l:139,54))=0.0 
;angled  portions  of  bow  and  stem 
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for  n=y0+l,y0+(9.0/23.0)*ny  do  begin 
templ=x0+(9.0/94.0)*nx-deltal 
temp2=x0+(89.0/94.0)*nx+delta2/3.0 
for  m=templ,temp2  do  begin 

if  outline(m,n)  le  BF*bkgdave(n)  then  outline(m,n)=0.0 
endfor 

boat(x0+(9.0/94.0)*nx-deltal,n)=255 

boat(x0+(89.0/94.0)*nx+delta2/3.0,n)=255 

outline(0:x0+(9.0/94.0)*nx-deltal-l,n)=0.0 

outline(x0+(89.0/94.0)*nx+delta2/3.0+l:139,n)=0.0 

deltal=deltal+l 

delta2=delta2+l 

endfor 

;top  of  angled  portion  of  bow 

for  n=y0+(9.0/23.0)*ny+l,y0+(15.0/23.0)*ny  do  begin 
templ=x0+(75.0/94.0)*nx 
temp2=x0+(89.0/94.0)*nx+delta2/3.0 
temp3=x0+(22.0/94.0)*nx 
temp4=x0+(58.0/94.0)*nx 
for  m=templ,temp2  do  begin 

if  outline(m,n)  le  BF*bkgdave(n)  then  outline(ni,n)=0.0 
endfor 

for  m=tenip3,temp4  do  begin 

if  outline(m,n)  le  BF*bkgdave(n)  then  outline(m,n)=0.0 
endfor 

boat(x0+(89.0/94.0)*nx+delta2/3.0,n)=255 
boat(temp3  ,n)=255 
boat(temp4,n)=255 
boat(templ,n)=255 

outline(x0+(89.0/94.0)*nx+delta2/3.0+l:139,n)=0.0 
outline(x0+(58.0/94.0)*nx+l:x0+(75.0/94.0)*nx-l,n)=0.0 
outline(0:x0+(22.0/94.0)*nx-l,n)=0.0 
delta2=delta2+ 1 
endfor 

boat(x0+(75.0/94.0)*nx:x0+(89.0/94.0)*nx+delta2/3.0,y0+(15.0/23.0)*ny)=255 
;fantail  and  foredeck 
temp5=5d)+(9.0/23.0)*ny 
boat(x0:x0+(22.0/94.0)*nx,temp5)=255 
boat(temp4:templ,temp5)=255 
;back  of  superstructure 
templ=x0+(22.0/94.0)*nx 
temp2=x0+(58.0/94.0)*nx 
temp4=y0+(  15 .0/23.0)*ny+ 1 
temp5=y0+(19.0/23.0)*ny 
temp6=y0+(9.0/23.0)*ny 
for  n=temp4,temp5  do  begin 
for  m=templ,temp2  do  begin 

if  outline(m,n)  le  duck*bkgdave(n)  then  outline(m,n)=0.0 
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endfor 

endfor 

boat(tenipl,temp4:temp5)=255 
outline(0;x0+(22.0/94.0)*nx-l,temp4;139)=0.0 
outline(x0+(58.0/94.0)*nx+l:139,temp4:139)=0.0 
;top  of  superstructure  rear  of  stack 
templ=x0+(22.0/94.0)*nx 
temp2=x0+(58.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny 
boat(teinpl:x0+(31.0/94.0)*nx,temp3)=255 
;rear  of  stack 
teinpl=x0+(31.0/94.0)*nx 
temp2=x0+(44.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny+l 
temp4=3K)+ny-l 
for  n=teinp3,temp4  do  begin 
for  m=templ,temp2  do  begin 

if  outline(m,n)  le  duck*bkgdave(n)  then  outline(m,n)=0.0 
endfor 
endfor 

boat(templ,temp3:temp4)=255 
outline(0:templ-l,temp3+l:139)=0.0 
outline(temp2+l:139,teinp3+l:139)=0.0 
;top  of  stack 

boat(templ:temp2,temp4)=255 
outline(*,temp4+l:139)=0.0 
;front  of  stack:  stack  average  already  calculated 
templ=x0+(44.0/94.0)*nx 
temp3=y0+(19.0/23.0)*ny+l 
temp4=y0+ny-l 
boat(temp  l,temp3  :temp4)=255 

;top  of  superstructure  front  of  stack:  already  calculated 
boat(x0+(44.0/94.0)*nx:x0+(58.0/94.0)*nx,y0+(19.0/23.0)*ny)=255 
;front  of  superstructiire:  superstructiure  already  calculated 
boat(x0+(58.0/94.0)*nx,y0+(9.0/23.0)*ny:y0+(19.0/23.0)*ny)=255 


•  ^  Hi ^  ^ ^  ^  :1c  ^ ^  ^ ^  ^ ^ ^  4:  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

;reposition  outline  if  needed 

window, 20, xsize=140,ysize= 140, xpos=150,ypos=0 
tvscl,boat 

read,'reposition  outline?  yes=l  no=0',repos 
if  repos  eq  1  then  goto,redo 

J 
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;calcTilate  results  and  display 

if  image  eq  1  then  begin 
;save  image  and  average  arrays 
bkgdavel=bkgdave 
outline  l=outline 
imagel=boat 

;increment  and  go  to  next  image 
image=image+l 
goto,nextimg 
endif 

if  image  eq  2  then  begin 
;save  image  and  average  arrays 
outline2=outline 
for  n=yO,3K)+ny  do  begin 
shptotl(n)=total(outlinel(*,n)) 
shptot2(n)=total(outline2(*,n)) 
indexl=where(outlinel(*,n),coimt) 
shppixl(n)=count 
index2=where(outline2(*,n), count) 
shppix2(n)=coimt 
if  shppix2(n)  ne  0.0  then  begin 

shpavel(n)=shptotl(n)/shppix2(n)  ;use  horizontal  ship  area  to  find  average 
shpave2(n)=shptot2(n)/shppix2(n) 
endif  else  begin 
shpavel(n)=0.0 
shpave2(n)=0.0 
endelse 
endfor 

bkgdave2=bkgdave 

image2=boat 

;calculate  average  ship  and  bkgd  values  above  and  below  horizon 
maxl=max(where(shpavel  gt  bkgdavel)) 
minl=min(where(shpavel  gt  bkgdavel)) 
shpbell=total(shpavel(minl:yvl))/(yvl-minl+l) 
shpabvl=total(shpavel(yvl+l:maxl))/(maxl-(yvl+l)+l) 
bkgbell=total(bkgdavel(minl:yvl))/(3rvl-minl) 
bkgabvl=total(bkgdavel(yvl+i:maxi))/(maxl-(3wl+l)+l) 

;display  image  with  outline 
title3=titlel 

window, 3, xsize= 140, ysize=140,xpos= 150  ,ypos=0,title=title3 
tvscl,imagel 

;generate  and  display  outlined  image  with  excess  background 
outlinel=outline(x0-10:x0+nx+20-l,y0-10:yfi+ny+20-l) 
shpimgl=boat(x0-10:x0+nx+20-l,y0-10:y0+ny+20-l) 
title5=titlel 

window,5,xsize=nx+20,ysize=ny+20,xpos=300,ypos=0,title=title5 

tvscl,shpimgl 
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;display  ship  and  background  average  profiles 
title7='shipl' 

window, 7, xsize=280,ysize=280,xpos=450,ypos=0,title=title7 

plotjShpavel 

oplot,bkgdavel 

;calculate  and  display  contrast  profile 
cntrstl=(shpavel-bkgdavel)/bkgdavel 
title9='cntrstl' 

window, 9, xsize=280,ysize=280,xpos=0,ypos=300,title=title9 
plot,cntrstl(niinl:maxl) 

;calculate  contrast  above  and  below  horizon  for  first  image 
cntlbel=(shpbell-bkgbell)/bkgbell 
cntlabv=(shpabvl-bkgabvl)/bkgabvl 
;calculate  average  ship  and  bkgd  values  above  and  below  horizon 
min2=min(where(shpave2  gt  bkgdave2)) 
max2=max(where(shpave2  gt  bkgdave2)) 
shpbel2=total(shpave2(min2:3rv2))/(yv2-min2+l) 
shpabv2=total(shpave2(yv2+l:max2))/(max2-(yv2+l)+l) 
bkgbel2=total(bkgdave2(min2:yv2))/(3rv2-min2+l) 
bkgabv2=total(bkgdave2(yv2+l:max2))/(max2-(yv2+l)+l) 

;display  image  with  outhne 
title4=title2 

window,4,xsize=140,ysize=140,xpos=150,ypos=160,title=title4 

tvscl,unage2 

;generate  and  display  outlined  image  with  excess  background 
title6=title2  .1^ 

outline2=outline(x0-10:x0+nx+20-l,y0-10:y0+ny+20-l) 

shpimg2=boat(x0-10:x0+nx+20-l,y0-10:y0+ny+20-l) 

window,6,xsize=nx+20,ysize=ny+20,xpos=300,ypos=160,title=title6 

tvscl,outline2 

;generate  and  display  impolarized  outlined  image  with  excess  backgroimd 
shpimgup=shpimgl/2.0+shpimg2/2.0 
outlinup=outlinel/2.0+outline2/2.0 
titleO='unp' 

window,0,xsize=nx+20,ysize=ny+20,xpos=300,ypos=80,title=title0 

tvscl,shpimgup 

;display  ship  and  backgroimd  average  profiles 
title8='ship2' 

window,8,xsize=280,ysize=280,xpos=730,ypos=0,title=title8 

plot,shpave2 

oplot,bkgdave2 

;calculate  and  display  horizontal  contrast  profile 
cntrst2=(shpave2-bkgdave2)/bkgdave2 
titlel0='cntrst2' 

window,10,xsize=280,ysize=280,xpos=0,yi)os=300,title=titlel0 

plot,cntrst2(min2:max2) 

;calculate  and  display  unpolarized  contrast  profile 
shpaveup=shpavel/2.0+shpave2/2.0 
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bkgaveup=bkgdavel/2.0+bkgdave2/2.0 

cntrstup=(shpaveup-bkgaveup)/bkgaveup 

titlel2='cntrstup' 

window, 12, xsize=280,ysize=280,xpos=0,ypos=0,title=titlel2 
plot,cntrstup(yO:yO+ny-l) 

;calculate  contrast  above  and  below  horizon  for  second  image 
cnt2bel=(shpbel2-bkgbel2)/bkgbel2 
cnt2abv=(shpabv2-bkgabv2)/bkgabv2 
;calculate  contrast  above  and  below  horizon  for  unpolarized  image 
shpupbel=shpbell/2.0+shpbel2/2.0 
shpupabv=shpabvl/2.0+shpabv2/2.0 
bkgupbel=bkgbell/2.0+bkgbel2/2.0 
bkgupabv=bkgabvl/2.0+bkgabv2/2.0 
cntupbel=(shpupbel-bkgupbel)/bkgupbel 
cntupabv=(shpupabv-bkgupabv)/bkgupabv 
;calctdate  contrast  improvement  above  and  below  horizon  for  v  and  h 
impbell=cntlbel/cntupbel 
impabvl=cntlabv/cntupabv 
impbel2=cnt2bel/cntupbel 
impabv2=cnt2abv/cntupabv 
;print  contrast  values  above  and  below  horizon 
print, 'contrast  above  horizon' 
print,'  vert  hor  map' 

print,cntlabv,cnt2abv,cntupabv 
print," 

print,'contrast  below  horizon 
print,'  vert  hor  imp' 

print,cntlbel,cnt2bel,cntupbel 
print," 

print, 'contrast  improvement  above/below  horizon' 
print, 'vertical:  ',impabvl,impbell 

print,  'horizontal: '  ,impabv2,impbel2 
endif 

;endfor 
goto, loop 
problem:  close,  1 

print, 'I  had  a  problem  and  so  I  quit' 
end 


90 


B.  CONTIGUOUS  AREA  METHOD  PROGRAM 


;flirashp  finds  the  ship  in  a  picture  completely  automatic 
;cursor  is  asked  for  if  at  picture  boimdary 
;stores  output  in  chosen  file  name 
;uses  shipfind 

;modular  version  uses  readflir  and  shipfind 

iscale=2  ;scale  displayed  pictures  by  2x 

close, 1  ;close  in  case  of  crash 

on_ioerror, problem  ;if  error,  quit  program 


pathl=' 

filel=' 

read,'input  path  to  *.img  files' ,pathl 
if  pathl  ne  "  then  path=pathl 
print, 'path=  ',path 

read,'input  type  of  files  eg  *.img  or  names  to  search  for',filel 

if  filel  eq  "  then  filel='*.img 
files=findfile(path+'  \  '+filel,count=num) 
for  i=0,num-l  do  print,files(i),i 

.Mc*******si=****]Qop  point  to  start  program  over  after  successful  processing**** 

loop:  print,n_elements(files)-l,'  last  number' 
read,'mput  i  for  number  to  plot  or  -1  to  stop',i 
istart=i 

if  i  It  0  then  goto,  problem 
read,'input  stop  number' ,istop 

ififix((istop-istart+l)/2)*2  ne  fix(istop-istart+l))  then  begin 
print,istart,istop,'  must  use  even  number  of  points' 
goto,  loop 
endif 

read,'input  first  guess  difference  t  level,  eg  5  to  include  all  of  ship  ,tdif 
tdifsave=tdif  ;save  tdif  so  that  changes  will  not  affect  batch 


.***********************2oop  point****************************** 

contrast=dblarr((istop-istart+l)/2,5) 
for  u=istart,  istop,2  do  begin 
for  i=ii,ii+l  do  begin 

tdif=tdifsave  ;restore  the  tdif  if  changed  in  loop 

f=files(i) 


file=f 

readflir,f,yy,t,har,tt,head  ;modular  input  returns  binary,  temparature 


.»»»»*******=t:*********Hi‘gp1ay  data  in  windows  for  visual  inspection*********** 


!order=0 
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ct=0 

loadct,ct 

iflii  eq  i)  then  window, 0,xsize=280,ysize=280,xpos=0,ypos=0,title='0  binary  data'  $ 
else  window, 5, xsize=280,ysize=280,xpos=280,ypos=0,title='5  binary  unsealed  data' 

print,'UNScaled  plot  of  binary  data  in  O' 

;tvscl,congrid(yy,140*iscale,140*iscale,Interp=l) 

tv,congrid(yy,140*iscale,140*iscale,Interp=l) 

wait,l 

loadct,16 

window,l,xsize=280,ysize=330,xpos=600,ypos=0,title='l  temperature' 
tvscl,rebin(tt,  140*iscale,  165*iscale) 


goto,skip 

print,'Cbck  in  scale  bar  to  exit' 

print, 'click  in  picture  for  point  temperature' 

replot:  cursor,xv,3rv,3/device 

if  yv  It  50  then  goto,skip 

print,t(xv,330-3rv),mint,maxt 

goto,  replot 

skip: 


.H:**^:**H:***M:****=i:*H:******|.gjj.g  baj-g  right  and  left  of  pictuTe  and 

^^6  loWBr  aVerOge  Value,  theU  SmOOth  ft 

thres0=(t(137,*)+t(138,*)+t(139,*))/3.0 

thres=(t(0,*)+t(l,*)+t(2,*))/3.0 

thres=thres  >  thresO  ;changed  to  larger  of  backgrounds 

minthres=min(thres) 

maxthres=max(thres) 

avthres=total(thres)/n_elements(thres) 

thres=smooth(thres,3) 

thres=rebin(thres,140,140)  ;make  the  vertical  line  xy  array 

thresO=thres  ;save  in  thresl  for  debug  140x140  array 


;ship=(t  gt  thres+tdif)*t  ;ship  is  greater  than  threshold 

;shipb=(t  gt  thres+tdif)*5^  ;get  binary  version  for  display 
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ship=(t  gt  tlires)*t 


;ship  is  greater  than  threshold 


shipb=(t  gt  thres)*yy  ;get  binary  version  for  display 

;new  version  uses  threshold  tdif  later 

ifltotal(ship)  eq  0)  then  stop  ;zero  values  mean  error-* 


;***make  linear  array  for  averages 

shipt=ship(where(ship  gt  0))  ;Pick  only  nonzero  values  for  average 

;shipt  is  linear  array  now 
avship=total(shipt)/n_elements(shipt) 
maxship=max(shipt) 
minship=min(shipt) 

;print,ininship,maxship,avship,'  ship  min, max, average  temperature' 
;rint,minthres,maxthres,avthres,  '  min  max  average  backgroimd' 


;Find  the  MAXIMUM  temperature 

point********************************************** 

maxt=max(ship,it)  ;find  the  max  in  the  ship  thresholded 


ix=it  mod  140 
iy=it/140 

window, 4,xsize=280,ysize=280,xpos=600,title=’#4  binary  ship’ 
tvscl,rebin(shipb,140*iscale,140*iscale)  ;make  ship  twice  as  big  for  display 

tvcrs,ix*iscale,iy*iscaleydevice  ;put  the  cursor  on  the  peak  location 

;print,'right  button,  use  auto  position,  left,  new  cursor  position' 

!err=4 

ifti  ne  istart)  then  begin 

print,ix,iy,nx,ny,xO,yO,format='(6(i3,2x),A19)','  ix,iy,nx,ny,x0,y0' 

iftix  It  10  or  iy  It  10)  then  cursor, x0,5K)/device  ;AUTOMATIC  LINE 

OUT 

iftnx/ny  gt  5  or  ny  It  10  or  nx  It  10)  then  cursor,x0,y0ydevice 
;AUTOMATIC  LINE  OUT  ;left  button  1,  middle  2  right  4 
endif 

iftlerr  eq  1)  then  begin 
iy=yO/iscale 
ix=x0/iscale 
x0=ix 
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yO=iy 


endif  else  begin 
xO=ix 


yO=iy 

endelse 


^  ^  1  1  .f^  1  *  j  *1^  xl> vXji  ^  ^  ^  ^  nl#  ^  ^  ^  ^  ^  ^  ^  ^  ■!»  mjm  ^L» xL* 

»*f;«j»^ jfj jf;5|5 jjjjt' hlT^  I  Tl n /  i  T^Om  r ’‘**^^ ^^*t..j.*i.«T.»f..x.«i.^«i..T..7.«r..i..f.*i»*i».t..f..t..x*'T"r*7*’T''i**T**i»*i»»T»v^'T*^ 


;following  necessary  for  reship  loop,  otherwise  ship  will  not  change  with  threshold 
reship: 

•  ^  H- ^  ^  ^  ^  ^  ^  ^  H*  H*  ^  ^  Jp j 


xO=ix 

yO=iy 


shipfind,shipb,ship,xO,yO,nx,ny,iscale  ;use  find  on  the  binary  version  to  avoid  neg  t 


;WAENING  SHIPFIND  SETS  ALL  VALUES  IN  SHIP  AND  SHIPB  ZERO  NOT  ON 
THE  SHIP 

tvscl,rebin(shipb,140*iscale,140*iscale)  ;show  the  contiguous  thresholded=0  ship 


;box:xy=boxfun(xO,yO,nx,ny/init,/message);option  to  readjust  the  box 

boxdraw,xO,yO,nx,ny  ;show  where  the  box  is 

xO=xO/iscale  ;rescale  from  picture  to  data  size 

3K)=y<)/iscale 

nx=nx/iscale 

ny=ny/iscale 


;print,x0,3d),nx,ny 

thres=(t(xO+nx,*)) 
thresl=t(xO,*) 
thres=thres  >  thresl 


;right  side  of  box 
;left  side  of  box 

;changed  to  greater  than  for  background 


thres=sniooth(thres,2) 

thres=rebin(thres,140,140)  ;make  an  array  140x140 

thresl=fltaiT(140,140) 

thresl(*,*)=0.0 

thresl(xO:xO+nx,yO:yO+ny)=thres(xO:xO+nx,yO:5K)+ny)  ;extract  part  of  array 
;Now  use  the  threshold  to  clean  up  the  ship  and  check  for  loop 
;the  ship  and  shipb  cannot  be  any  larger  than  tdif=0,  so  whittle 
ship=(t  gt  thresl+tdiD*ship  ;ship  is  greater  than  threshold 

shipb=(t  gt  thresl+tdiD*shipb 
;get  binary  version  for  display 
x0=ix 


94 


yO=iy 

shipfind,shipb,ship,xO,yO,nx,ny,l  ;use  find  on  the  binary  version  to  avoid  neg  t 


window, 2, xsize=280,ysize=280,xpos=300,title='#2  local  threshold' 
tvscl,rebin(shipb,140*iscale,140*iscale)  ;redisplay 

;above  is  entire  140x140  picture  thresholded 

box=t(x0:x0+nx,y0:y0+ny)  ;temperature  in  box 

raw=yy(x0:x0+nx,y0:y0+ny)  ;raw  signal  in  box 

sMpl=ship(xO:xO+nx,yO:yO+ny)  ;ship  temperature  in  box 

rawl=shipb(x0:x0+nx,y0;y0+ny)  ;ship  binary  in  box 

shiplb=  box*(shipl  eq  0)  ;ship  temperatme  BACKGROUND  in  box 

rawlb=  raw*(rawl  eq  0)  ;ship  binary  BACKGROUND  in  box 
;shipl  and  rawl  are  box  sized  subset  of  ship  and  shipb 

window, 3, xsize=nx*iscale*2,ysize=ny*iscale*2,xpos=0,ypos=0 


tvscl,congrid(shipl,(nx+l)*iscale*2,(ny+l)*iscale*2) 

tdifl=0.0 

;read,'input  new  tdif  or  0  to  stop  ',tdifl  ;AUTOMATIC  LINE  OUT 
ifitdifl  gt  0)  then  begin 
tdif=tdifl 

;remake  the  ship  variables  using  the  NEW  box  threshold 

ship=(t  gt  thres+tdif)*t  ;ship  is  greater  than  threshold 

shipb=(t  gt  thres+tdif)*yy  ;get  binary  version  for  display 


goto,reship 


;loop  up  to  reship 


endif 

.:|:********************JJggJJJp  LOOP 

point************************************************* 

stuff***************************************** 

ifii  eq  ii)  then  begin 


endif 


tdifsave=tdif 

x0save=x0 

y0save=y0 

nxsave=nx 

nysave=ny 

ixsave=ix 

iysave=iy 


;box  position 


;peak  position 
;peak  position 


if(i  eq  ii)  then  begin 
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shiplsv=shipl 

rawlsv=rawl 

shiplbsv=shiplb 

rawlbsv=rawlb 

tsave=t 

yysave=yy 

sbipsv=ship  ;save  tbe  temperature  tbresbolded  full  pictxire 

sbipbsv=sbipb  ;save  tbe  binary  tbresbolded 

endif  else  begin 

;if  tbe  first  box  is  bigger,  renormalize  to  it 

;otherwiser  renormabze  all  data  to  the  second  box 

;The  following  are  necessary  for  correct  contrast  using  binary 

;rawl,rawlb,rawlsv,rawlbsv  -  only  two  of  these  are  destroyed 

;only  one  set  of  the  two  needs  to  be  recomputed  based  on  the  new  box 

;size  which  is  larger 

print,  nxsave,nysave,nx,ny,'  nxsave,y,nx,ny’ 
ifinxsave*nysave  gt  nx*ny)  then  begin 
print,'using  old  ship  size' 

;use  save  box  and  make  rawl  and  rawlb  match  old  ship  shipsv 

xO=xOsave 

yO=yOsave 

nx=nxsave 

ny=nysave 

rawl=yy*(shipbsv  gt  0)  ;mask  saved  correct  ship  on  cmrent  data 

rawlb=yy*(shipbsv  eq  0)  ;mask  for  background  =  not  ship 

rawl=rawl(x0:x0+nx,y0:y0+ny)  ;select  the  current  (saved  )  box 
rawlb=rawlb(x0:x0+nx,y0:y0+ny) 


endif  else  begin 

print, 'using  new  ship  size' 

rawlsv=  yysave*(shipb  gt  0)  ;mask  current  ship  on  old  data  140x140 
rawlbsv=yysave*(shipb  eq  0) 

;current  not  ship  binary  times  old  data 

rawlsv=rawlsv(x0:x0+nx,y0:y0+ny)  ;now  select  the  current  box 
rawlbsv=rawlbsv(x0:x0+nx,y0:y0+ny)  ;backgrovmd 
;it  would  also  be  nice  later  to  update  the  shipl  (temperature)  stuff 
endelse 
endelse 


endfor 
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> 

loadctjO 
;;in  order 

;ship  1  average  value 

;ship  1  background  average  value 

;ship  2  average  value 

;sliip  2  background  average  value 

;ship  unpolarized 

;ship  background  unpolarized 

;ship  unpolarized  average  value 

;ship  background  luipolaxized  average  value 

;contrast  shipl 

;contrast  shipsv  (shipO) 

;contrast  unpolarized 

;contrast  shipl/contrast  unpolarized 

;contrast  shipsv/contrast  unpolarized 

shiplav=total(shipl)/total((shipl  gt  0)) 
shiplbav  =  total(shiplb)/total((shiplb  gt  0)) 
shipsav=total(shiplsv)/total(  (shiplsv  gt  0)) 
shipsbav  =  total(shiplbsv)/total(  (shiplbsv  gt  0)) 
shipup=(shipl+shiplsv)/2.0 
shipbup=(shiplb+shiplbsv)/2.0 

shipupav=total(shipup)*2.0/total((shipl  gt  0)+(shiplsv  gt  0)) 
shipubav=total(sliipbup)*2.0/total((shiplb  gt  0)+(shiplsv  gt  0)) 

;above  are  taken  for  temperature,  below  for  raw  dataq 

shiplav=total(rawl)/total((rawl  gt  0)) 
shiplbav  =  total(rawlb)/total((rawlb  gt  0)) 
shipsav=total(rawlsv)/total(  (rawlsv  gt  0)) 
shipsbav  =  total(rawlbsv)/total(  (rawlbsv  gt  0)) 

;WAENING,  adding  binary  will  exceed  255  and  give  error 
;one  must  first  float  the  niunbers;  using  temperature  above  is  ok 
shipup=(float(rawl)+float(rawlsv))/2.0 
shipbup=(float(rawlb)+float(rawlbsv))/2.0 
shipupav=total(shipup)*2.0/total((rawl  gt  0)+(rawlsv  gt  0)) 
shipubav=total(shipbup)*2.0/total((rawlb  gt  0)+(rawlbsv  gt  0)) 


;total  shipup  gt  0  does  not  work  because  of  registration 


cshipl=(shiplav-shiplbav)/shiplbav 
cshipsv=(shipsav-shipsbav)/shipsbav 
cshipu  =(shipupav-shipubav)/shipubav 
cclu=cshipl/cshipu 


ccsu=cshipsv/cshipu 

;rawl  =  ship  picture  in  binary  0-255  (in  box) 

;rawlsv=  same  picture  0 

;rawlb  =  binary  backgroimd  in  box 

;rawlbsv  =  binary  saved  backgroimd 

;yy  =  full  picture  (same  as  t  for  temperature 

;yysave=  picture  0  binary 

;t=temperature  pictmre 

;tsave=old  picture  temperature 

J 

Print, 'Contrast  and  enhancement, first  file  =',file,tdif,tdifsave 
ifiistop-istart  gt  0)  then  begin 

print, files(ii+l),files(ii),'  unpolarized  ,first/impolarized,second/impolarized' 

print,cshipl,cshipsv,cshipu,cclu,ccsu 

contrast((ii-istart)/2,0:4)=[cshipl,cshipsv,cshipu,cclu,ccsu] 

;save  the  values 

endif  else  begin 

;print,'  first  contrast  ',file 

;print,cshipl 

endelse 


jdisplay  difference  image  for  both  polarizations 
window,2,xsize=nx*4,ysize=ny*4,xpos=0,3ipos=0,title=files(istop) 

tvscl,rebin(rawl,(nx+l)*4,(ny+l)*4) 

window,3,xsize=nx*4,ysize=ny*4,xpos=0,5^os=140,title=files(istart) 


tvscl,congrid(rawlsv,(nx+l)*4,(ny+l)*4) 
^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 


in  automatic 


goto,  noplot 

;plot  ship  average  value  compared  to  background  average  value  and  make  .tif  file 


;horizontal  polarization 
titlel='horizontal'+'  polarization' 

window,4,xsize=350,ysize=300,xpos=300,3T}os=0,title=titlel  ;files(istop) 


wset,4 

;call  to  external  function  avrow  averages  rows 
plot,avrow(shiplsv),indgen(140),background=255,color=0 

oplot,avrow(shiplbsv),indgen(140),color=150 
;  oplot,bkgdave(*),indgen(140),color=150 
plotl=tvrd() 
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plotl=plotl(*,299-indgen(300)) 
tiff_write, title  1+'  .tif  ,plotl 

;plot  comparison  of  ship  average  values  with  respect  to  polarization  and  make  .tif  file 
title2='polarization'+'  comparison' 

window, 5, xsize=350,ysize=300,xpos=300,ypos=330,title=title2  ;'vertical/horizontar 
wset,5 

plot,avrow(shipl),indgen(140),background=255,color=0 

oplot,avrow(shiplsv),indgen(140),color=150 

plot2=tvrd() 

plot2=plot2(*,299-indgen(300)) 

tifF_write,title2+'.tif,plot2 

;plot  first  ship  average  value  compared  to  background  average  value  and  make  .tif  file 

;vertical  polarization 
title3='verticar+'  polarization' 

window, 6, xsize=350,ysize=300,xpos=660,3T>os=0,title=title3  ;files(istart) 


wset,6 

plot,avrow(shipl),indgen(140),background=255,color=0 

oplot,avrow(shiplb),indgen(140),color=150 

plot3=tvrd() 

plot3=plot3(*,299-indgen(300)) 

tifF_write,title3+'.tif,plot3 

;plot  comparison  of  backgrounds  with  respect  to  polarization  and  make  .tif  file 
title4='background'+'  comparison' 

window,7,xsize=350,ysize=300,xpos=660,3^os=330,title=title4 


wset,7 

plot,avrow(shiplbsv),indgen(140),backgrormd=255,color=0 


oplot,avrow(shiplb),indgen(140),color=150 

plot4=tvrd() 

plot4=plot4(*,299-indgen(300)) 
tiff_write,title4+'  .tif  ,plot4 
.*H:H:********H:************gjj^(j  gf  temperatuTe 


windows************************************* 


;plot  ship  average  value  compared  to  backgroimd  average  value  and  make  .tif  file 
;horizontal  polarization 
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titlel='hzpor 

window, 10, xsize=350,ysize=300,xpos=300,ypos=0,title=titlel  ;files(istop) 

;call  to  external  function  avrow  averages  rows 

plot,avrow(rawlsv),indgen(140),background=255,color=0,title='horizontar+'  polarization' 

oplot,avrow(rawlbsv),indgen(140),color=150 
;  oplot,bkgdave(*),indgen(  140),color=150 
plotl=tvrd() 

plotl=plotl(*,299-indgen(300)) 

tiff_write,titlel+'.tif,plotl 

;plot  comparison  of  ship  average  values  with  respect  to  polarization  and  make  .tif  file 
title2='comp' 

window,12,xsize=350,ysize=300,xpos=300,57pos=330,title=title2 

;'vertical/horizontar 

plot,avrow(rawl),indgen(140),backgroimd=255,color=0,title='polarization'+' 

comparison' 

oplot,avrow(rawlsv),indgen(140),color=150 

plot2=tvrd() 

plot2=plot2(*,299-indgen(300)) 

tiff_write,title2+'.tif,plot2 

;plot  first  ship  average  value  compared  to  background  average  value  and  make  .tif  file 

;vertical  polarization 
title3='vpor 

window,13,xsize=350,ysize=300,xpos=660,ypos=0,title=title3  ;files(istart) 

plot,avrow(rawl),indgen(140),background=255,color=0,title='verticar+'  polarization' 

oplot,avrow(rawlb),indgen(140),color=150 

plot3=tvrd() 

plot3=plot3(*,299-indgen(300)) 

tifiF_write,title3+'.tif,plot3 

;plot  comparison  of  backgrounds  with  respect  to  polarization  and  make  .tif  file 
title4='bkgnds' 

window,14,xsize=350,ysize=300,xpos=660,ypos=330,title=title4 

plot,avrow(rawlbsv),indgen(140),background=255,color=0,title='backgroimd'+' 

comparison' 

oplot,avrow(rawlb),indgen(140),color=150 

plot4=tvrd() 

plot4=plot4(*,299-ir  ?n(300)) 
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tiff_write,title4+'.tif,plot4 


noplot: 

endfor 

;saveflr8.pro  prints  results 
;saveflr8.pro  prints  results 
fileout='  ' 

close, 2 

read, 'input  output  filename',fileout 
openw,2,fileout 

printf,2,systime(),'  date  of  analysis' 
printf,2,'Temperature  difference  =  ',tdif 
printf,2,istop-istart+l,'nuniber  of  files' 
j=istart 

for  i=0,(n_elements(contrast(*,0))-l)  do  begin;printf,2,files(2*i),files(2*i+l) 
printf,2,files(j+l),files(j),fonnat='(2(A30))’ 

j=j+2 

printf,2,'first,second,unpolarized  ,first/unpolarized,second/unpolarized' 

printf,2,contrast(i,*),format="(5(f9.6,5x))" 

print,i,contrast(i,*),fonnat="(i3,lx5(f9.6,5x))" 

endfor 

close,2 

print, 'finished' 

goto,  loop 
problem:  close,  1 

print, 'I  had  a  problem  and  so  I  quit' 
end 


C.  BOX  METHOD  PROGRAM 

;lasthope.pro 

close,  1 
intfii=0 

shpnmbr=:l 

loadct,ct 


!order=0 

print,'flirs  makes  movie  of  Airs' 

;close  in  case  of  crash 

ct=0 


;wmdow,0,xsize=280,ysize=280,xpos=0,3q)os=0 
on_ioerror, problem 


pathl=' 

filel=' 


'  prmt,'file  size  must  be  20446 

with  header'  read,'input  path  to  *.img  files' ,pathl 

if  pathl  ne  "  then  path=pathl 
print, 'path=  ',path 


101 


f 


read, 'input  type  of  files  eg  *.img  or  names  to  search  for',filel 
if  filel  eq  "  then  filel='*.img 

files=findfile(path+'  \  '+filel,count=num)  for 

i=0,num-l  do  print, files(i),i 

^  *1#  ^ ^  ^ ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

•  5ft  5f5  5Jt  JJt  5p  ^  5|t  5^  5|^  ^  ^  5^  ^  5ft  ^  ^  ^  "i^  ^  ^  ^  ^  ^  ^  ^  ^  ^  #y% ^  ^ ^  ^ #j» #f»  rj»  »f%  pjm  0j^  •f*  ^  vjwvj* 

J 

;print  time  of  analysis  and  number  of  files  analyzed 
> 

printf,7,'  ',systime() 

printf,7,'  ',strmid(files(5),17,5),'.txt' 


loop:  print,n_elements(files)-l,'  last  number' 

read,'input  i  for  number  to  plot  or  -1  to  stop',i  if  i  It  0 

then  goto,  problem  istart=i 

read, 'input  stop 

munber',istop  for  i=istart,  istop  do  begin 

f=files(i) 


9 

;open  output  file 

9 

openr,l,f 

print,'opening  ',f,i 

head=bytarr(846) 

y=bytarr(140,70) 

yl=bytarr(140,70) 

yy=bytarr( 140, 140) 

readu,l,head 


close,  1 
index=indgen(70) 
yy(*,index*2)=yl(*,index) 
yy(*,index*2+ l)=y(*, index) 
ifimax(yy)  eq  255)  then  begin 
yl=yy(  where(yy  gt  254)) 
nover=n_elements(yl) 


endelse 

yy=yy(*,139-indgen(140)) 


readu,l,y 

readu,l,yl 


endif  else  begin 
nover=0 


9 

;save  ship  images 

shipl=b5d;arr(140,140) 
ship2=bytarr(140, 140) 
if  shpnmbr  eq  1  then  begin 
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shipl=yy 
shipll=yy 
shpnmbr=2 
shpname  l=files(i) 
title  l='vert' 

window, l,xpos=0,3^os=0,xsize= 140, ysize=140,title=titlel 
tvscl,shipl 
endif  else  begin 
ship2=yy 
shpname2=files(i) 
title2='hor' 

window,2,xpos=0,ypos=160,xsize=140,ysize=140,title=title2 

tvscl,ship2 

shpnmbr=l 

endelse 

endfor 

y= 


yy=congrid(37y,140*2,140*2,Interp=l) 

loadct,0 

;wset,0 

print,'Scaled  plot  of  data  in  O' 

range=double(head,770) 

level=double(head,762) 

dlevel=double(head,618) 

eniiss=double(head,626) 

ct=double(head,650) 

r=double(head,570) 

b=double(head,578) 

f=double(head,586) 

aint=double(head,674) 

at=double(head,666) 

pl=emiss*ct 

p2=(  l-emiss)/emiss 

p3=1.0-ct 

iamb=r/(exp(b/anit)-f) 

iatm=r/(exp(b/at)-f) 

y=range*(float(yy)-128.)/256.0+level+dlevel 
y=y/pl-p2*iamb-p3/pl*iatm 
t=b/alog(abs((r/y+f)))-273 . 15 
idate=head(482:487) 
itime=head(488:495) 

;  get  rid  of  old  variables 
;delvar,y,yl,yy,head 


yl=l 

;!order=l 

;tvscl,yy 


103 


shpup=fltarr(  140) 

shpup(*)=0.0 

shpl=fltarr(  140) 

shpl(*)=0.0 

shp2=fltaiT(  140) 

shp2(*)=0.0 

bkgdup=fltarr(  140) 

bkgdup(*)=0.0 

bkgdl=fltarr(140) 

bkgdl(*)=0.0 

bkgd2=fltarr(  140) 

bkgd2(*)=0.0 

shpave=fltarr(  140) 

shpave(*)=0.0 

;shpavel=fltarr(140) 

;shpavel(*)=0.0 
;shpave2=fltarr(  140) 

;shpave2(*)=0.0 
;shptotl=fltarr(  140) 

;shptotl(*)=0.0 
;shptot2=fltarr(  140) 

;shptot2(*)=0.0 

shppixl=fltarr(  140) 

shppixl(*)=0.0 

shppix2=fltarr(  140) 

shppix2(*)=0.0 

outline=fltarr(140,140) 

outline(*,*)=0.0 

outline  l=fltarr( 140, 140) 

outlinel(*,*)=0.0 

outline2=fltaiT( 140, 140) 

outline2(*,*)=0.0 

image=l 

nextimg: 

wset,image 

zz=boxfun(xO,yO,nx,ny,init=intfn/message) 

intfh=l 

box=tvTd(x0,3K),nx,ny) 
bkgdave=fltarr(  140) 
bkgdave(*)=0.0 
shppix=fltaiT(  140) 
shppix(*)=0.0 


;assign  image  array  to  outline  array  and  mark  horizon 
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J 


if  image  eq  1  then  begin 
outline=shipll 

print, 'click  on  horizon  or  bottom  to  exit' 
replotl:  cursor, xv,yv,3ydevice 
if  yv  It  5  then  goto,exitl 
xvl=xv 
yvl=yv 
yvy=yv 

print, 'horizon=',5rvl 
print, 'click  on  point  or  bottom  to  exit' 
goto,replotl 
exitl:print,'exiting' 
endif 

if  image  eq  2  then  begin 
outline=ship2 

print,'click  on  horizon  or  bottom  to  exit' 
replot2:  cursor,xv,3^,3,/device 
if  yv  It  5  then  goto,exit2 
xv2=xv 
yv2=yv 
yvy=yv 

print, 'horizon=',yv2 
print,'click  on  point  or  bottom  to  exit' 
goto,replot2 
exit2:print,'exiting' 
endif 

5 

;calculate  backgrovmd  average  values 

J 

if  5K)+ny/2.0  It  69  then  begin 
;box  below  narcissus 
;background  average  below  box 
for  n=0,3K)-l  do  begin 
bkgdave(n)=total(outline(0:139,n))/140.0 
endfor 

;backgroimd  average  aroimd  box 
for  n=3K),3K)+ny-l  do  begin 

bkgdave(n)=(total(outline(x0-5:x0-l,n))+total(outline(x0+nx:x0+nx+4)))/10.0 

;background  average  above  box 
for  n=3K)+ny,139  do  begin 

bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 
endif  else  begin 


endfor 


endfor 
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;box  above  narcissus 
;background  average  below  box 
for  n=0,y0-l  do  begin 

bkgdave(n)=(total(outline(0:19,n))+total(outline(120:139,n)))/40.0 

;background  average  around  box 
for  n=yO,yO+ny-l  do  begin 

bkgdave(n)=(total(outline(x0-5:x0-l,n))+total(outline(x0+nx:x0+nx+4)))/10.0 

;background  average  above  box 
for  n=y0+ny,139  do  begin 
bkgdave(n)=(total(outline(0:139,n)))/140.0 
endfor 
endelse 

read, 'enter  background  factor', dd 
duck=dd 

if  image  eq  1  then  begin 
factorl=duck 
endif  else  begin 
factor2=duck 
endelse 

again: 

shppix(*)=0 

outline(0:139,0:y0-l)=0.0 
outline(0:139,y0+ny:139)=0.0 
outline(0:x0- 1,0 : 139)=0 . 0 
outline(xO+nx:  139,0: 139)=0.0 
for  n=y0,y0+ny-l  do  begin 
for  m=x0,x0+nx-l  do  begin 

if  outline(m,n)  le  duck*bkgdave(n)  then  outline(m,n)=0.0 
if  outline(m,n)  gt  0.0  then  shppix(n)=shppix(n)+l 
endfor 

if  shppix(n)  eq  0  then  shpave(n)=0.0 

if  shppix(n)  ne  0  then  shpave(n)=tot8d(outline(*,n))/shppix(n) 
endfor 

numbel=total(shppix(yO:yvy)) 
numabv=total(shppix(5wy+ 1 :54)+ny-l)) 

shpabv=total(outline(x0:x0+nx-l,yvy+l:3K)+ny-l)) 

shpbel=total(outline(xO:xO+nx-l,yO:yvy)) 

bkgdabv=total(bkgdave(yvy+l:yO+ny-l))/((yO+ny-l)-(yvy+l)+l) 
bkgdbel=total(bkgdave(yO  :yvy  ))/(yvy-y0+ 1) 

window, 3,xsize=140,ysize=140,xpos=420,ypos=0 
tvscl,outline 


endfor 


endfor 
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read, 'enter  new  background  factor  or  zero  to  continue', qq 
duck=qq 

if  duck  gt  0.0  then  begin 
if  image  eq  1  then  begin 
outline=shipll 
factorl=duck 
goto,again 
endif  else  begin 
outline=ship2 
factor2=duck 
goto,again 
endelse 
endif 

if  image  eq  1  then  begin 
outlinel=outline 
shpl(*)=shpave(*) 
bkgdl(*)=bkgdave(*) 
shpabvl=shpabv 
shpbell=shpbel 
bkgdabvl=bkgdabv 
bkgdbell=bkgdbel 
numabvl=numabv 
numbel  l=mmibel 
image=image+ 1 

window,4,xsize=140,ysize=140,xpos=280,5q)os=0 

outline(*,3wl)=255 

tvscl,outline 

window,  5  ,xsize=  140,ysize= 140  ,xpos= 140,ypos=0 
tv,shipll 
goto,nextimg 
endif 

if  image  eq  2  then  begin 
outline2=outline 
shp2(*)=shpave(*) 
bkgd2(*)=bkgdave(*) 
bkgdup(*)=bkgdl(*)/2.0+bkgd2(*)/2.0 
shpup(*)=shpl(*)/2.0+shp2(*)/2.0 

shpabv2=shpabv 

shpbel2=shpbel 

bkgdabv2=bkgdabv 

bkgdbel2=bkgdbel 

shpabvla=shpabvl/numabvl 

shpbella=shpbell/numbell 

numabv2=numabv 

numbel2=numbel 
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if  numabv2  ge  ntnnabvl  then  begin 
shpabvl=(shpabvl+bkgdabvl*(nmnabv2-nuniabvl))/nuniabv2 
shpabv2=shpabv2/nuniabv2 
endif  else  begin 
shpabvl=shpabvl/numabvl 

shpabv2=(shpabv2+bkgdabv2*(numabvl-nuniabv2))/numabvl 

endelse 

if  nuinbel2  ge  numbell  then  begin 
shpbell=(shpbell+bkgdbell*(nunibel2-nunibell))/nunibel2 
shpbel2=shpbel2/nunibel2 
endif  else  begin 
shpbell=shpbell/numbell 

shpbel2=(shpbel2+bkgdbel2*(nuinbell-nmnbel2))/nunibell 

endelse 

cntabvl=(shpabvl-bkgdabvl)/bkgdabvl 

cntbell=(shpbell-bkgdbell)/bkgdbell 

cntabv2=(shpabv2-bkgdabv2)/bkgdabv2 

cntbel2=(shpbel2-bkgdbel2)/bkgdbel2 

shpabvu=shpabvl/2.0+shpabv2/2.0 

shpbelu=shpbell/2.0+shpbel2/2.0 

bkgdabvu=bkgdabvl/2.0+bkgdabv2/2.0 

bkgdbelu=bkgdbell/2.0+bkgdbel2/2.0 

cntabvup=(shpabvu-bkgdabvu)/bkgdabvu 

cntbelup=(shpbelu-bkgdbelu)/bkgdbelu 

impabvh=cntabv2/cntabvup 

inipbelh=cntbel2/cntbelup 

inipabvv=cntabvl/cntabvup 

inipbelv=cntbell/cntbelup 


window, 6, xsize=140,ysize=  140, xpos=280,ypos=160 

outline(*,yv2)=255 

tvscl,outline 

window,7,xsize=140,ysize=140,xpos=140,3T)os=160 

tv,ship2 

endif 


print,'  shpbella  shpabvla 

print,shpbella,shpabvla 

print,'  shpbell  shpbel2  shpbelu' 

print,shpbell,shpbel2,shpbelu 

print,'  bkgdbell  bkgdbel2  bkgdbelu' 

print,bkgdbell,bkgdbel2,bkgdbelu 
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print," 

print,'  cntabvl  cntabv2  cntabvup' 

print,cntabvl,cntabv2,cntabvup 

print,'  cntbell  cntbel2  cntbelup' 

print,cntbell,cntbel2,cntbelup 

print,'  numabvl  nuniabv2' 

print,nuniabvl,nuniabv2 

print,'  numbell  nunibel2' 

print, nxxmbel  1  ,nunibel2 

print,'  impabw  impabvh' 

print,impabw,impabvh 

print,'  impbelv  impbelh' 

print,impbelv,impbelh 

print,'bkgd  factor  v/h',factorl,factor2 

print,'horizon  v/h',yvl,5rv2 

ship=string('ship') 

bkgd=string('bkgd') 

cntabv=string('cntabv') 

cntbel=string('cntber) 

numabv=string('nuniabv') 

numbel=string('niunber) 

impabv=string('impabv') 

impbelsstringCimpbel') 

BF=string('BF') 

hrzn=string('hrzn') 

printf,?,'  ',shpnamel,'\',strmid(shpname2,17,12) 

printf,7,'vert','hor','nnpor,itime(0),itime(2),itime(4),format='(20x,a4,5x,a3,6x,a5,2x,i2,lh:, 

i2,lh:,i2)' 

printf,7,'shipabv',shpabvl,shpabv2,shpabvu,format='(10x,a7,lx,fB.2,3x,fB.2,3x,fB.2)' 

printf,7,'bkgdabv',bkgdabvl,bkgdabv2,bkgdabvu,format='(10x,a7,lx,fB.2,3x,fB.2,3x,fB.2)' 

printf,7,'shipber,shpbell,shpbel2,shpbelu,format='{10x,a7,lx,fB.2,3x,fB.2,3x,fB.2)' 

printf,7,'bkgdber,bkgdbell,bkgdbel2,bkgdbelu,format='(10x,a7,lx,fB.2,3x,f6.2,3x,fB.2)' 

printf,7,'cntabv',cntabvl,cntabv2,cntabvup,format='(10x,a6,2x,fB.2,3x,fB.2,3x,fB.2)' 

prmtf,7,'cntber,cntbell,cntbel2,cntbelup,format='(10x,a6,2x,fB.2,3x,fB.2,3x,fB.2)' 

printf,7,'nuniabv',numabvl,numabv2,format='(10x,a6,2x,fB.2,3x,fB.2)' 

printf,7,'number,nunibell,numbel2,format='(10x,a6,2x,fB.2,3x,fB.2)' 

printf,7,'impabv',impabw,impabvh,format='(10x,a6,2x,fB.2,3x,fB.2)' 

printf,7,'impber,impbelv,inipbelh,format='(10x,a6,2x,fB.2,3x,fB.2)' 

printf,7,'BF',factorl,factor2,fonnat='(10x,a2,6x,fB.2,3x,fB.2)' 

printf,7,'hrzn',yvl,3^2,format='(10x,a4,4x,fB.2,3x,fB.2)' 

printf,7," 


intfh=0 

shpnnibr=l 
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goto,  loop  problem: 

close,  1  print,'!  had  a 

problem  and  so  I  quit' 
end 
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